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EARTHQUAKE RESPONSE OF FIXED OFFSHORE STRUCTURES 


Owen Major Kirkley, Ph.D. 
Department of Civil Engineering 
University of Illinois at Urbana-Champaign, 1973 

The response of fixed offshore structures, such as those used for 
oi] drilling and production purposes, to earthquake excitation is investigated. 
An approximate method of computing maximum response quantities in the structure 
is developed. The approximate method uses superposition of classical normal 
mode responses as determined by response spectrum techniques. The structure 
is modeled as a lumped parameter system, the earthquake is represented by the 
deterministic time history record of the ground acceleration, and the hydro- 
dynamic forces are considered to be those computed from Morison's equation. 

A single degree of freedom oscillator with nonlinear drag character- 
istics 1s eee to various representations of ground motion, including 
earthquake motions. The existence of a high frequency tail in the response 
Spectrum is shown, and a method for linearizing the nonlinear hydrodynamic 
drag force is developed and shown to be satisfactory. A method is developed 
for conservatively estimating the magnitude of equivalent linear hydraulic 
damping in the system based on velocity characteristics of the ground motion. 
Rules are developed for the construction of idealized response spectra for 
several different response quantities. The effects of the nonlinear hydro- 
dynamic drag force on the inelastic response of the single degree of freedom 
oscillator are investigated. 

The principles and approximate methods developed during the study 


of the single degree of freedom system are applied to the multi-degree of 





freedom model of the actual structure. The linearization technique is ex- 
tended to the multi-degree of freedom system in such a manner as to insure 
the existence of classical normal modes in the equations of motion. A Series 
of multi-degree of freedom systems is subjected to one of the earthquake 
records and the linearization technique is shown to be satisfactory in com- 
puting the response quantities of interest. The response of a series of 
offshore towers to an earthquake record is computed by the superposition of 
modal responses as determined by both exact spectra and idealized spectra. 
The approximate method considers the effects of the structure-fluid inter- 
action by conservatively estimating the amount of linear hydraulic damping 


developed as a result of this interaction. 





CHAPTER | 
INTRODUCTION 
1.1 General 


The design of fixed offshore structures is frequently based on the 
static forces resulting from the loads imposed on the structure by wind, 
waves, ocean currents, and earthquakes. Design on the basis of static forces 
is satisfactory when the response of the structure is not magnified by the 
dynamic nature of the loads. However, the characteristics of the structure 
and the forces acting upon it are often such that the response of the struc- 
ture to wind, waves, currents, and earthquakes is dynamic and magnified greatly 
over the static response. In these cases, design should properly be based on 
the dynamic response of the structure. Various methods are available for cal- 
culating the dynamic response of a land-based structure to earthquake forces, 
but similar calculations for offshore structures are complicated by the pres- 
ence of nonlinear hydrodynamic forces acting on submerged portions of the 
structure. Further, these nonlinear forces are dependent not only on the 
physical characteristics of the structure, but also on the actual response of 
the structure to the forces acting on it. 

Many offshore structures are currently in use, with the most notable 
use being that of oi] drilling and production platforms. While platforms in 
water depths of up to 300 feet are fairly common now, the search for new 
sources of oi] and gas may result in many offshore structures situated in 
greater water depths. The location of oil resources off the coast of California, 


in the Gulf of Alaska, and other portions of the circumpacific earthquake belt 





as well as other regions of seismic activity emphasizes the need to be able 

to predict the response of offshore structures to earthquake motions. Drill- 
ing and production platforms represent a very large investment in cost, and 
large numbers of personnel are involved in the operation of such a structure. 
Failure of one of these platforms could result in many lives being lost as 
well as the loss of investment both from the cost of the platform and the loss 


of production time. 
1.2 Description of the Problem 


An offshore structure of the type mentioned in the preceding sec- 
tion is shown schematically in Fig. 1.1(a) and is a continuous system with 
an infinite number of degrees of freedom. Such structures are space frames 
composed of circular pipe sections generally welded together, with diagonal 
bracing as shown in the figure. The top portion of the structure may consist 
of more than one level and contains the heavy equipment necessary for 07] 
drilling and production and often contains facilities for berthing and messing 
of the crew. The operating and production portion of the structure is gener- 
ally quite heavy. The research reported in this thesis is based on a lumped 
parameter model!” of the structure with mass, stiffness, forces and other 
characteristics of the structure lumped into a discrete number of nodes, as 
shown in Fig. 1.1(b). 

The equations of motion for the system shown in Fig. 1.1(b), sub- 


jected to ground motion in quiescent water, are, in matrix form: 


PmJix.} + (C.Jtx, - x} + [K]ix,-x,} = {P} Gap 


Xx 
g 


Numbers refer to List of References at end of text. 





where 


mass matrix of the system 
structural damping matrix of the system 
stiffness matrix of the system 


absolute acceleration of node j 


absolute velocity of node j 


absolute displacement of node j 


ground velocity 


ground displacement 


The symbol [ -~] denotes a matrix in which the off-diagonal elements are zero. 


No structural dampers are shown on the model, but the damping matrix fea can 


take a variety of forms as will be discussed in a later chapter. 


The hydrodynamic forces acting at each node of the system are con- 


sidered to be those computed from Morison's equation? for circular cylinders 


as: 


where 


mix} - Pky Tk |x, |} (1.2) 


added mass matrix? = e[~C VJ 
mass density of the fluid surrounding the structure 
empirical inertia coefficient at node j 


enclosed volume at node j 


drag factor matrix = sey A. 





cD 
i 


D empirical drag coefficient at node j 


> 
i 


projected area at node j 


If the coordinate transformation, {Z} = {X. - Kats 1s made, where 
{Z} is a vector of a motion quantity of a node relative to that of the ground, 


Eqs. 1.1] and 1.2 can be combined to yield: 


' 


CM J{Z} + [C1{Z} + [K]}{Z} = -PMJ{%} - espe led 


where 


virtual (or total) mass matrix = Pmj+fm 1 


bud 


‘ 
g 


ground acceleration 


Equation 1.3 differs from the equations of motion for a land-based 
structure subjected to earthquake forces in two important aspects. First, 
the presence of the added mass term causes the structure to possess natural 
frequencies and mode shapes which are different from those of a structure with 
the same stiffness and only the structural mass [~m.J. Secondly, the hydro- 
dynamic force term Kp ~H{x. |X,|} is not present in the conventional land-based 
earthquake problem, and this nonlinear hydrodynamic force causes the major dif- 
ficulty in solving Eq. 1.3. 

If the nonlinear hydrodynamic force [Kp tx, |x,|} can be converted 
into an “equivalent" linear hydrodynamic. force, say Be palate EQ. |. 35¢an 


be written as 


DM J{Z} + (C12) + CKI{Z} = -DM IK} - Pky I} (1.4) 


Since {Z} = {X, ~ Xqis Eq. 1.4 can be rewritten as 





DMZ} + [C.1(Z} + Ok TZ} + CKI{Z) 
= - DMC} - [kp JOxg) a5) 


Equation 1.5 differs from the equations of motions for a land-based multi- 
degree of freedom damped structure of mass [>M.], subjected to a base exci- 
tation, in the linear drag force T Kp tx, }- 

One well-known and convenient method of solving simultaneous linear 
differential equations of motion is the normal mode method.’ In the normal 
mode method, the vector of a response quantity expressed in one coordinate 
system, such as the relative displacement between each node and the base {Z}, 


1s transformed into the "normal" coordinate system by the relation 


7A Lo ){Y) (1.6) 
where 
[¢] = a matrix of the natural mode shapes of the system with the 
ith column in the matrix representing the deflected shape 
of the ith mode of the system. 
and {Y} = a vector of modal response quantities with the ith element 


in the vector representing the modal response quantity in 


the ith mode. 


Thus, Eq. 1.6 expresses the response at each node in the model as a combina- 
tion of the responses at that node in each of the natural modes of vibration 

in the system. The natural mode shapes in the modal matrix [9] are orthogonal?. 
with respect to [-M.] and [K] and will be orthogonal with [C.] if the elements 
of [C.] are specified in a certain way.*! Substitution of Eq. 1.6 into Eq. 1.5 


yields, 





CMJCo](¥} + CC JLo]{¥} + Op ~ILo]{v} + CKILo{Y} 
= - PMC} - Pp Six) (1.7) 
Pre-multiplying Eq. 1.7 by [4]! yields, 
Co} CMILolt¥} + Co]'Cc.Iont¥} + [o]'C Ky, ~ILolt¥} 
+ Lo] TEKILO}3 = -Lo]E MIG} - Lely TF (1.8) 


Recognizing the orthogonal relationships between the modal matrix [o] and 


other matrices in Eq. 1.8, the following matrices are defined: 
[-M*_] = Co) RM Jo] = a diagonal modal mass matrix 


me =] - (ol'Cc. IC] = a diagonal modal structural damping 
matrix. The form of [C.] which permits 


this is discussed later. 


CK J = Co) KIC] 


AY 


[Kp] = Col DK Jo] 


a diagonal modal stiffness matrix 


a modal linear hydraulic damping matrix 


which is not necessarily diagonal 
= * 
we] - mM A TARY = a matrix of circular natural frequencies 


of the system 
Using these quantites, Eq. 1.8 becomes 
ee ee fe ° n, 2 * 
CM JY} + DC ~}{Y} + [Kp HY} + [DK YY} 


= -Le]'D MIG} - [o]'E Ky IK) (1.9) 





If the modal linear hydraulic damping matrix [Ky] can be made dia- 
gonal, then the set of simultaneous differential equations of motion repre- 
sented by Eq. 1.9 will be completely uncoupled. Suppose that the nondiagonal 
matrix [Ky] can, in some satisfactory manner, be converted into a diagonal 


matrix, say _ Then Eq. 1.9 can be written 
* ee * Pe * @ * 
eee ec iy + TK, Jit + Pk ty) 
= -[6]'D MIG} - [e]'DKp Tig) (1.10) 


Or, for the ith mode, 


fee C OUYS eoOk, Y. + KY 
ii ij Sa, i DL. . j lade" 
" a . 
= = o.. M..1(X_) - te (xe) era ae) 
j=. Ji JJ} 9 j=l DL, g 
where 
nn = number of nodes in the model 
Equation 1.11 can also be written as 
ere ene. + Qe. Ge Yo ao. Y 
7] Ss. 1 7j h. ji 1 177 
nn nn 
+ ial 2, 93% So, 
= =% Jee e 
= . z (Xq) = ‘ (XQ) (1.12) 
Ma 11 
where Cc" 
Be = aoe = structural damping ratio in the ith mode 
7 M. .w 
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wis 
; * 
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K 
= linear hydraulic damping ratio in the ith mode 


Bh 
If the modal response quantity of interest, Y, Y, or Y can be ob- 
tained from Eq. 1.12 for each mode, then the total response of the structure 
can be obtained by using the relationship given in Eq. 1.6. The modal re- 
sponse quantities vary with time, and the time variation is different for 
each mode, so an exact calculation of the total response of the structure 
can only be made if the exact time histories of the modal responses are known 
and superimposed on one another. Generally, the maximum response quantity 
at each node of the model is desired, but the modal maxima do not necessarily 
occur at the same time. An approximation to the actual maximum response at 
each node of the model can be made by determining the maximum response in each 
mode and then combining the modal responses in some fashion. The response 


spectrum is a convenient way to determine the maximum modal response.?‘~}? 


1.3 Summary of Previous Research 


Equation 1.3, or a similar form, has been solved by other investi- 
gators uSing various methods. Much of the existing work in this area has 
been devoted to the response of offshore structures to wave forces, but the 
nonlinear difficulties are somewhat similar for wave forces and earthquake 
forces acting on submerged or partially submerged structures. Further, the 
bulk of recent investigations has been devoted to the use of stochastic repre- 
sentations of wave and earthquake motions. 

Burke and Tighe” integrate the equations of motion directly for a 


series of towers subjected to different deterministic ground acceleration 





records by using an iterative technique which computes the response quanti- 
ties of interest at step-by-step time points throughout the earthquake record. 
For a tower in 800 feet of water subjected to four accelerograms (normalized 
to the intensity of the 1940 El Centro earthquake), the effect of the water- 
structure interaction was found to reduce the maximum base shear and the maxi- 
mum overturning moment by about 33 percent from the values found when the 
water-structure interaction was not included in the solution. This method, 
while permitting the direct use of the structural properties of the system, 
including the nonlinear drag force, requires the use of a digital computer 
and the availability of actual or simulated earthquake records. For purposes 
of preliminary design of a structure for a "design" earthquake of the designer's 
choice, however, the step-by-step technique is costly, time consuming, and 
requires that a time history of the "design" earthquake be available. Pre- 
liminary design would best be served by a rapid, approximate method of solu- 
tion which permits the use of selected parameters for the design earthquake 
and accounts for the effects of the nonlinear drag forces, yet yields satis- 
factory results. 

Penzien, in association with other investigators,° ° linearizes 
the drag forces in Eq. 1.3 by replacing the nonlinear force vector with a 
linear force vector which performs the same amount of work during the history 
of the response. This linearization results in an "equivalent" linear hydraulic 
damping matrix acting on the system in much the same manner as the structural 
damping matrix acts. The linear equations of motion are then transformed into 
a normal coordinate system and solved by the normal mode method of superposition.’ 


Unfortunately, the linear hydraulic damping matrix does not uncouple when the 
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equations are transformed into the normal coordinate system. As a result, 

an uncoupling (or diagonalizing) technique, similar to the original linear- 
izing technique, must be employed in order to separate the equations into the 
desired modal equations. The linearizing and diagonalizing are accomplished 
by an iterative process which uses actual response quantities of the system 
as they are developed, and the procedure is time consuming and requires great 
care in its formulation and use. For towers subjected to a stochastic repre- 
sentation of the ground motion, equivalent linear hydraulic damping ratios on 
the order of 9 percent of critical are reported for the first two modes of 
towers in deep water, with "very much smaller" hydraulic damping ratios for 
the tower in shallow water. ° 

Foster? solves the problem of an offshore tower subjected to 
stochastic representation of wave forces by linearizing the drag force, then, 
recognizing that nonproportional hydraulic damping may result, solving the 
linear equations of motion using complex eigenvalues and eigenvectors. This 
procedure recognizes that, unless the linearization of the drag force is done 
in a certain way, classical normal modes may not necessarily exist in the equa- 
tions, thus the normal mode method of superposition is not a rigorously cor- 
rect approach. 

Nath and Harleman!® conclude that, for large diameter cylinders in 
deep water waves, the nonlinear drag force can be neglected, and, since shallow 
water structures are relatively stiff, and therefore the dynamic response to 
waves is not greatly magnified over the static response, design of shallow water 
Structures can proceed on a static basis. These conclusions, however, were 
arrived at by considering the response of structures to waves, which generally 


have much longer period components of motion than do earthquakes. 
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1.4 Objective and Scope 


The objective of this thesis is to investigate the effects of the 
nonlinear hydrodynamic drag force on the earthquake response of single degree 
of freedom and multi-degree of freedom systems and to develop a satisfactory 
method of linearizing this force in order that conventional methods of solving 
linear differential equations may be employed in the solution of the equations 
of motion. 

A further objective is to develop an approximate method for the 
solution of Eq. 1.3 by the superposition of normal mode responses, using 
response spectra to estimate the modal response quantities necessary. A rapid, 
yet fairly accurate method of solution is desirable in the preliminary stages 
of the design of an offshore structure in order that the designer may have 
an estimate of the response of a tentative design to one or more design earth- 
quakes of his specification. It should be kept in mind that the approximate 
solution developed in this thesis is no substitute for a more detailed analysis 
of the structure once the design is in its final stages. 

In Chapter 2, the effects of the nonlinear hydrodynamic drag force 
on the response of a single degree of freedom oscillator are explored. A 
reasonable method of linearizing the drag force and a method for conservatively 
predicting the resultant linear drag effects are developed and shown. A method 
is developed for estimating certain response quantities of the oscillator in 
order to determine more accurately the range of hydraulic drag and inertia 
coefficients which should be used in the force equation. The response of the 


oscillator in the inelastic region, assuming elasto-plastic force-displacement 


er 
a 





characteristics of the spring, is shown. Finally, the construction of ideal- 
ized response spectra, which include the effects of hydraulic drag, is shown. 
The concepts developed in Chapter 2 for the single degree of freedom 
system are extended to the multi-degree of freedom system in Chapter 3. The 
nonlinear drag force is linearized in such a manner as to insure the existence 
of classical normal modes in the equations of motion in order that normal mode 
superposition may be used as an approximate solution to the problem. The method 
of linearization is shown to be accurate for a series of five degree of free- 
dom systems which are representative of the type of structure being investi- 
gated. A means of estimating the amount of equivalent linear hydraulic damping, 
using ground motion velocity characteristics, is developed and is shown to 
yield good results. The construction of an idealized response spectrum for 
each mode of the multi-degree of freedom system is discussed and described, 
and sample problems are solved by the superposition of normal mode responses. 
The application of the approximate method to a series of actual off- 
Shore towers is demonstrated, and the results of calculations of the response 


of these towers are discussed in Chapter 4. 
1.5 Related Topics 


In this section, several topics of general interest in the field of 
earthquake engineering are briefly discussed in order that the reader may more 
fully understand the display of results by the use of response spectra, the 
numerical method of integrating the equations of motion, and the characteris- 


tics of the earthquake records used in the research. 
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ieoet Response Spectra 


The equation of motion for a damped, single degree of freedom os- 


cillator in air, subjected to a base disturbance is 

m+ C(x, - Xq) + K(x, - x9) = 0 (1.13) 
By making the transformation Z = ee Xq> Eq. I.13 becomes 

mZ + C2 + KZ = - m Te 
Or, in an alternate form 


Z + 28. wl + we Z = pe (1.15) 


The solution of Eq. 1.15 depends on the characteristics of the ground motion 
function, the amount of damping in the system as represented by the dimension- 
less parameter Be» and the natural frequency of the system, w. The frequency 
Parameter enters into the solution as a dimensionless product of w and a 
Characteristic time of the base motion. ?*? 

A plot of the maximum value of a response quantity, without regard 
to sign, as computed from Eqs. 1.13, 1.14 or 1.15, as a function of the fre- 
quency of the system is called a response spectrum for that quantity. The 
response spectrum is a convenient way to represent various maximum response 
quantities of interest as a function of the variables involved in the solution. 
The response quantity of greatest interest is usually the maximum spring dis- 
tortion, PA nis which occurs during the response history of the oscillator, 
although other response quantities are sometimes of interest. ?? 


The response spectrum is conveniently shown on tripartite logarithmic 
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graph paper as shown in Fig. 1.2. The grid of this representation is con- 
structed such that the abscissa represents the frequency of the oscillator, 
the ordinate represents velocity units of response, the diagonal lines 
sloping up to the right are lines of constant displacement units, and the 
lines sloping down to the right are lines of constant acceleration units. 
In this thesis, unless otherwise indicated, the response spectra shown are 


deformation response spectra in which the pseudo-velocity, w|Z| is plotted 


max 
against the undamped frequency of the system. The displacement quantity 


shown is the maximum spring deformation, |Z|___, and the acceleration quantity 


max 


is the pseudo-acceleration, w°|Z| Response spectra for a damped single 


max 
degree of freedom oscillator, in air, are shown in Fig. 1.2 for the El Centro 
earthquake. The shape of the spectra is typical for strong motion earthquakes 
of the type used in this thesis. More detailed information on response spectra 
can be found in the literature.*?~** 

Since the response spectra for strong motion earthquakes exhibit the 
characteristic shape shown in Fig. 1.2, the spectra can be idealized by straight 
line segments. Amplification factors, which relate the idealized spectrum to 
have been developed 


the ground motion quantities x » and xX 


X 
g max’ “g max g max’ 


for a range of damping ratios.'* The ability to construct an idealized repre- 
sentation of the response spectrum for a given set of conditions permits an 
estimate of the response to be made without having an actual earthquake record 


available. 


1.5.2 Numerical Integration of the Equations of Motion 


The equations of motion for a system subjected to a base disturbance 





Is 


can be solved by numerical integration. In this thesis, numerical integra- 
tion of the equations of motion is performed by the Newmark "Beta Method."® 

In the calculations, 8 is set to be equal to 1/6, thus representing the ground 
response acceleration as linear between adjacent time points. The time inter- 
val used is 1/20 of the natural period for frequencies up to 10 cps and 1/40 
of the natural period for frequencies equal to, or greater than 10 cps. These 
time increments insure good accuracy, stability and convergence in the solu-' 
tion. Convergence of the solution at any time point is assumed to have oc- 
curred if two successive iterations yield values of acceleration which are 
within 0.1 percent. In the case of multi-degree of freedom equations of 
motion, the time increment is based on the shortest natural period in the 
system, and iteration is carried out until the computed acceleration at each 
node is within 0.1 percent of the value computed on the previous cycle. In 
all cases, the integration of the equations of motion is carried out for a 
period of time equal to one-half the natural period (the longest natural 
period in the case of multi-degree of freedom systems) after the base distur- ~ 


bance ends. 


1.5.3 Earthquake Records Considered 


Four actual earthquake records are used in this thesis. The records 
are described numerically in Table 1.1 and the time history of the ground ac- 
celeration, velocity, and displacement for each record is shown in Figs. 1.3 
through 1.6. In each earthquake record, a very small (less than .005g) arti- 
ficial acceleration point has been generated near the end of the record so 


that the velocity of the base motion is zero at the end of the record. 





16 


Although a parabolic baseline correction’® had previously been applied to the 
acceleration records, a slight residual velocity still remained at the end of 
the record and it is this residual velocity which is removed by the artificial 
acceleration point. The four records chosen are considered to be representa- 


tive of strong motion earthquakes on firm soil. 
1.6 Notation 


As much as possible, the same symbols are used for both the single 
degree of freedom oscillator and the multi-degree of freedom system, with 
matrix and vector notation indicating the multi-degree of freedom system. 
Symbols are defined as they first appear in the text. For ease of reference 


the symbols used are summarized in this section. 


[~A.] projected area matrix 
Ch empirical dimensionless drag coefficient 
C, empirical dimensionless inertia coefficient 
[C.J structural damping coefficient matrix 
‘a : — dj 
mie modal damping coefficient matrix = [@] [C_Ifo] 
D diameter of a single degree of freedom oscillator 
{HDPF } vector of participation factors for hydrodynamic 


parameter response quantities 


{HOPF} = [K 97! Coq! cKIC7} 


{I} unit vector 





MJ 
cM 
{P} 


{PF} 


VJ 


{Y}, {Y}, {Y} 
{Qi}, {Q3, {0} 


mr, {7}, {72} 


vs 


structure stiffness matrix 
modal stiffness matrix = Co) CKI[0] 
drag factor matrix = se DCA 


modal drag factor vector for very high frequency 
ae 
system = DM J 'Co]'D Ket} 


linear drag factor matrix 


modal linear drag factor matrix = [ol'T Ky, ~IL0] 


dimensionless drag parameter 

total structure dynamic mass matrix = [mJ + Cm —) 
modal mass matrix = fo] RM JLo] 

hydrodynamic force vector 


participation factor vector for displacement response 


_ ell 
SN A aed) el She 
enclosed volume matrix 


displacement, velocity , and acceleration vectors 
in normal coordinates 


vectors of displacement, velocity, and acceleration 
of the structure relative to the base 


deformation at which yielding occurs 





Cm] 
tm ~] 
nm 


nn 


{xg} {xghs {xg} 


x X X 
g max’ “g max’ xg max 


{x}; {x}; {X.} 
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relative acceleration between fluid and a single 
degree of freedom oscillator 


Subscripts indicating mode number 
subscripts indicating node number 
structural mass matrix 
added mass matrix e[~C,V~] 

number of modes in the multi-degree of freedom system 


number of nodes in the lumped parameter model 


hydrodynamic force acting on the single degree of 
freedom oscillator 


relative velocity between fluid and a single degree 
of freedom oscillator 


ground displacement, velocity and acceleration vectors 


maximum displacement, velocity, and acceleration in 
the ground motion record 


vectors of structure displacement, velocity, and 
acceleration 


displacement, velocity, and acceleration of the ith 
mode single degree of freedom equation of motion 
after the modal participation factor has been taken 
Out. Expressed in normal coordinates 





{y} 


[¢] 


ig 


constant of proportionality between the linear drag 
factor matrix and the total mass matrix [Ky ~ 
= af MJ 


linear hydraulic damping ratio in the ith mode 


structura] damping ratio in the ith mode 
CPF i: 
g _ _s§ ii 


S e M 
] 2 e i). e 


vector of modal stiffness factors {fy} = Co) CK]{1} 
ductility factor for inelastic response 
max 
Z 
y 
mode shape matrix 


mass density of the fluid 


circular natural frequency of the ith mode 


Poe = eMeter kK] 
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CHAPTER 2 
SINGLE DEGREE OF FREEDOM OSCILLATOR 
2.1 General 


The normal mode method of solving Eq. 1.3 requires the solution 
of single degree of freedom equations of motion. For this reason, and to 
understand more fully the effects of the nonlinear hydrodynamic drag force 
on the response of a simple system, Chapter 2 is devoted to the investigation 
of the effects of the hydrodynamic drag force on the response of a single de- 
gree of freedom oscillator. 

In the investigation of the response of the oscillator with non- 
linear drag characteristics it is shown that the response of the oscillator 
in the very high frequency region can be estimated if the drag properties of 
the system are known and if the maximum velocity of the base motion is known. 
The nonlinear drag force is linearized by equating the work done by a ficti- 
tious linear force with the work done by the real, nonlinear force. This 
method of linearization is shown to be accurate until the high frequency re- 
gion of the spectrum is reached. It is shown that the linear drag force is 
equivalent to a linear hydraulic damper in which the force is dependent upon 
the velocity of the mass, and that the magnitude of the hydraulic damping 
coefficient is dependent upon the amount of structural damping in the system. 
A means of estimating the amount of linear hydraulic damping is developed 
for earthquake motions, and the estimate is shown to be conservative except 
in the low frequency range. The maximum response of the system in the mid- 


frequency range to earthquake excitation is shown to be dependent upon the 
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total amount of linear damping present, and not upon the ratio in which the 
damping is distributed between structural and hydraulic. 

Response quantities which may be used in order to select reasonable 
values of the inertia and drag coefficients, Cy and Cy» are explored, and a 
method of estimating these response quantities is developed. The inelastic 
response of an oscillator with elasto-plastic force-displacement character- 
istics of the spring is investigated and the results are shown in the form of 
response spectra. 

The construction of idealized elastic response spectra which include 
the effects of the nonlinear hydrodynamic forces is demonstrated in this 


chapter. 
2.2 Elastic Deformation Response Spectra 


The equation of motion can be written for the single degree of 


freedom oscillator as 


<e 














MZ + CZ + KZ = “MK - KX. ; (221) 
or 
aa ae ee (2.2) 
Sue ‘is g M “s!"'s ; 
or 
zs : 2 KD x max i ane 
Z + 2Bwl + wZ = -% - (— SM) (-e) «|x| 
S g K X s' Ss 
g max 
[23 ) 


As pointed out in Chapter 1, the solution of the equation of motion for a 


land-based single degree of freedom oscillator (Eq. 2.1 without the nonlinear 
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hydrodynamic drag force KnXe Xe |) depends on the characteristics of the ground 
motion, the amount of damping in the system as represented by the dimensionless 
parameter Bo» and the natural frequency of the system, w. The solution of 
Eq. 2.1 depends on the magnitude of Ky in addition to these parameters. A 


dimensional analysis of Eq. 2.1 yields the dimensionless drag parameter 


KD *g max 


K 


as one of the parameters on which the response of the oscillator depends. 
When using a maximum ground acceleration of 0.5 g, values of 0.01, 0.2, and 
0.6 represent a reasonable range over which to vary the dimensionless drag 
parameter. For the maximum ground acceleration in the EI] Centro earthquake 
record (0.312 g), values of 0.006, 0.08, and 0.50 for the parameter yield 
Slightly greater variations in the magnitude of Kp: Obviously, a zero value 
of the drag parameter represents a "no-drag" system and increasing values of 
the parameter represent increasing magnitudes of Ky in the system for fixed 
values of K and Xq ee 
Equation 2.1] is solved numerically for the steady-state maximum 


response of the oscillator when the base acceleration is 


Xg = 16.1 sin 2nt 


where 
Xo is in ft/sec“ 


t is in seconds 


Deformation response spectra (in which 2 | a is plotted against frequency) 


are plotted in Figs. 2.1 and 2.2 to show the results of the calculations. 
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In order to investigate the response of the oscillator to a slightly 
more complicated representation of the base disturbance, the base acceleration 
is represented as a six term cosine series. The numerical characteristics of 
the series used in the response calculations are shown in Table 2.1. The 
series has the features of zero acceleration, velocity, and displacement at 
the beginning and end of the base motion. Equation 2.1 is solved numerically 
for the maximum response of the oscillator when subjected to five cycles of 
this base disturbance, and the results are shown in the form of deformation 
response spectra in Figs. 2.3 and 2.4. 

Equation 2.1 1s solved numerically for the maximum response of the 
oscillator when subjected to the El Centro earthquake record. The results 
are shown in the form of deformation response spectra in Figs. 2.5 through 
e7. 

A study of the spectra presented in Figs. 2.1 through 2.7 reveals 
several interesting characteristics which are best discussed in terms of 
frequency ranges of the oscillator. 

In the low frequency region, the maximum spring deformation, |Z| 


max’ 


approaches that for the case of KX ay 


hydrodynamic drag in the system for the range of drag parameter used. This 


/K = 0, regardless of the amount of 


1s reasonable to expect, since the magnitude of the drag force is a function 
of the oscillator velocity, xe and very soft (low frequency) systems exhibit 
little velocity response of the mass (x.). 

In the mid-frequency region, where spectra for the conventional 
“no-drag" oscillator show an amplification of the pseudo-velocity over the 
maximum ground velocity, the spectra for systems with drag also show an ampli- 


fication of the pseudo-velocity over the maximum ground velocity. However, 
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the presence of the nonlinear hydrodynamic drag force causes a reduction in 
the response of the oscillator from the response of the "no-drag" oscillator. 
The amount of reduction is greater for greater values of the drag parameter 

up to some frequency. This reduction would suggest that, in the mid-frequency 
range of the spectrum, the nonlinear hydrodynamic force can be replaced by 

an equivalent linear hydrodynamic force which will effect the response in 

much the same manner as the presence of structural damping effects the re- 
sponse. The amount of reduction in response caused by the drag force is also 
seen from the spectra to be dependent upon the amount of structural damping 

in the system. This is easily noted by comparing Figs. 2.5 and 2.6. In Fig. 
2.5, the system has no structural damping, and a slight amount of drag (repre- 


sented by K.x /K = 0.006) is seen to reduce significantly the maximum 


D*g max 
response of the oscillator in the frequency range of about |] to 10 cps. In 
Fig 2.6 the system has 2 percent structural damping, and the drag parameter 

of 0.006 is seen to cause very little reduction in the response spectrum be- 
tween | and 10 cps. The effects of structural damping on the amount of re- 
duction in response caused by the hydrodynamic drag are to be expected, since 
the drag force is a function of the oscillator velocity, Xe and this velocity 
is a function of the amount of structural damping in the system. 

In the high frequency region, response spectra for "no-drag" oscil- 
lators typically show the pseudo-acceleration approaching the maximum ground 
acceleration, Xq _—_ Response spectra in the high frequency region for the 
oscillators with nonlinear hydrodynamic force exhibit entirely different 


Characteristics, however. In all the response spectra shown, a frequency is 


reached at which the maximum response quantity ee for the oscillator with 
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drag exceeds that for the "no-drag" oscillator. In fact, for the systems 
with larger values of the drag parameter, the quantity ele. is seen to ap- 
proach a relatively constant value. The magnitude of this high frequency 
"tail" in the spectrum can be estimated as follows: 


Consider an undamped system with spring stiffness so great that 


The forces on the oscillator are 


F spring i. F inertia i Farag 

pring - Ee u Kp%5 [Xs | 
or, 

os NOt Ke Yo 

spring g Dg 


Since the drag parameter 


K. X 
Dg max 
K 


is held constant, Ky must increase as the spring stiffness K increases, so 
a frequency will be reached at which the inertia force can be neglected with 


respect to the drag force. If the inertia force is neglected, then 


| ; 
een nax - Ky iXq nae 
And 
” °2 
Z| = F Spring! max : ‘0 (; Ne - (‘D *g max (9 MAX ) (2.4) 
max K K ‘“g max K X : 


g max 
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A comparison of [Zhe obtained at high frequencies by the numerical 
integration of Eq. 2.1 with (Zine computed from Eq. 2.4 is shown in Table 
2.2 for all systems and base disturbances considered. A study of the data 
presented in Table 2.2 shows that, particularly for the earthquake records, 
the high frequency tail of the deformation response spectrum can be satisfac- 
torily estimated by the use of Eq. 2.4. 

When an acceleration record is integrated numerically in order to 
obtain the velocity and displacement records, the integration is done from 
time point to time point in the record. This double integration yields the 
values of X age 12.79 in./sec and Xa es 8.24 inches for the El Centro 
record as shown in Table 1.1]. However, if the acceleration record is inte- 
grated to find the ground velocity and displacement at some time point which 
1s between time points in the acceleration record, linear interpolation of 
the acceleration record must be performed in order to obtain the acceleration 
quantity necessary to compute the ground velocity and displacement at the 
intermediate time point. This interpolation introduces small errors into 
the acceleration, which in turn introduces slightly larger errors into the 
ground velocity which in turn introduces even larger errors into the ground 
displacement computed at the intermediate time point.’® It was found that, 
when using the very small time increments necessary for stability and con- 
vergence in the numerical integration of the equation of motion at high fre- 
quencies (above about 8 cps), the interpolation of the El Centro acceleration 
record introduced errors into the magnitude of X, such that X max Was found 
to be about 14.0 in./sec in lieu of the value 12.79 in./sec previously ob- 


tained, even though double precision was used in all calculations. In order 
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to achieve comparability of results in the comparison of |Z|_., estimated from 


Eq. 2.4 and eae obtained at high frequencies from the integration of the 


equation of motion, . = 14.0 in./sec is used in Eq. 2.4 to obtain the 


*g max 


results shown in Table 2.2 for Kx /K = 0.006 and 0.08. 


D*g max 


2.3 Linearization of the Drag Force 


It was seen in Chapter 1 that the use of the normal mode method to 
solve Eq. 1.3 depends on the ability to linearize the nonlinear hydrodynamic 
force vector PK pix. |x, |} in a satisfactory manner. The linearization of 
the hydrodynamic drag force in the single degree of freedom system is inves- 
tigated in this section. It is desirable to replace the nonlinear drag 


DL*s? which will 


yield the same maximum response in the system. Consider the work done by the 


force, Kx. [Xe » in Eq. 2.1 by a linear drag force, say, K 


two forces:?° 


t 





Work done by the nonlinear _ 2 fe : 
drag force over a period of time | (Kx, 1X, | ) (x, dt 
O 
1b 
Work done by a linear drag 2 nity 
force over a period of time | (Kp) X) (x, dt 


Equating the amounts of work done, 


t 
| Kyx2 
O 


X 
> 








t °? 
pX< dt = | Koi %< dt 
0 


K <x |x Py = VK <xe> 
Se S 


D DL 


<xe |X |> 

SS 

ky | a | 
SX 67 


DL (2.5) 
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(The symbol < > denotes the average value of the quantity enclosed over the 
time of the response of the oscillator.) This expression is the same expres- 
sion used by other investigators®-? to linearize a nonlinear force. With the 


nonlinear drag force replaced by an "equivalent" linear drag force, Eq. 2.1 


becomes, 

MZ + C2 + KZ = “MX, - KoL%s 

MZ + CZ + Ky, Z + KZ = “MK - KDLXg (2.6) 
or, 

oe @ . 2 7 a ° 

Zot 28 wl + 28, wl + wZ = a 28 UX. (227) 
where, 

ADL 
By = on equivalent linear hydraulic damping ratio 


It now remains to investigate how the maximum response in the equi- 
valent linear system compares with the maximum response in the actual nonlinear 
system. 

Equation 2.6 is solved numerically for Pale using a linear drag 
factor Koy which is linearized in accordance with Eq. 2.5 for the single de- 
gree of freedom oscillator subjected to the series representation of the base 
disturbance. The quantities xe |X [> and <xo> necessary for the linearization 
were obtained from the integration of the nonlinear equation of motion for the 
same set of conditions. The results are plotted on Figs. 2.8 and 2.9, and it 
is seen that, for this base disturbance, the linearization is very good until 


the high frequency tail begins to appear. This might have been expected, 
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since a linearization of a weak nonlinearity would normally be better than a 
linearization of a strong nonlinearity, and the high frequency tail is the 
result of the nonlinear hydrodynamic drag force becoming quite strong and 
dominant in Eq. 2.1. 

A similar set of calculations is made for the El Centro earthquake. 
The response spectra for K, and K , are Shown in Fig. 2.10 for B. = 0.02 and 


D D 
Xx /K = 0.50. It is seen from this figure that the linearization for 


Dg max 
the earthquake motion, while not quite as good as the linearization for the 
series representation of the base motion, is still reasonably good up until 
the high frequency tail of the spectrum is reached. A comparison of Za 
computed for the nonlinear drag system with alee for the linear drag sys- 
tem iS shown in Figs. 2.11 through 2.13 for all of the variables considered. 
It is seen from a study of these figures that the maximum response of the 
linear system is close to the maximum response of the nonlinear system until 
the high frequency region is reached. The figures further indicate that the 
linearized response is slightly on the conservative side of the actual re- 
Sponse until the high frequency region is reached. The accuracy of the 
linearization in the high frequency region is expected to become poor as the 
nonlinearity is becoming more and more dominant in the system, and thus the 
linearization becomes less and less equivalent. However, Eq. 2.4 permits a 
fairly good estimate to be made of the maximum relative displacement in the 
high frequency region. 

Since the linearization of the drag force in Eq. 2.1 is seen to 
yield a slightly conservative maximum response at frequencies below the be- 


ginning of the high frequency tail in the response spectrum, a method of 
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estimating the linear drag factor KO from the known nonlinear drag factor 
Ky is now developed for earthquake motions. A study of Eq. 2.5 shows that 
two velocity quantities, <x |X |> and cxe> are necessary in order to linearize 
the drag force and, since these quantities are actual response quantities 
of the nonlinear system, they cannot be calculated a priori. A series of 
computations is made in which the quantity 

<x |x, [> 


2 
=X o> 
S 





is computed for the actual, nonlinear oscillator subjected to the El Centro 
earthquake. The results are shown in Figs. 2.14 through 2.16. From these 
figures it is seen that, for the range of variables studied, if the quantity 
0? e 
<xo |X. [> 


2 
<X > 
S 


is estimated as 


x |> 


“2 
<X > 
g 


Ze 





for the particular earthquake considered, then this estimated value will be 


less than the actual 





2? eS 
<x_ |X |> 
S S 
22 
“x > 
S 


for all frequencies except low frequencies. Using the value of 
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<x |x |> 
<> 
in Eq. 2.5 to linearize the drag term results in a value of KOE which is less 
than the value of Koy which would result from using the actual oscillator 
velocity quantities. When this lower value of Ko is converted into the 
linear hydraulic damping ratio, By» shown in Eq. 2.7, the value of By is 
less than the actual value of By which would be computed had the actual os- 


cillator velocity terms been available for use in Eq. 2.5. In other words, 


linearizing the drag term by using 


«xe |x |> 


<X_ > 
g 


in Eq. 2.5 leads to a conservative value of Kor which, in turn, leads to a 
conservative value of the maximum response in all except the low frequency 
region and the high frequency tail region of the spectrum. The unconserva- 
tive nature of Koy in the low frequency region when computed from the ground 
velocity terms is reasonable to expect, however, as the amount of hydraulic 
damping is dependent on the velocity response of the oscillator, and very 
soft (low frequency) systems have little velocity response, so the amount of 
hydraulic damping should approach zero as the system becomes softer and softer. 
Figures 2.14 through 2.16 further show that, for a system with given drag 


characteristics, the linear drag factor Kor estimated by using the quantity 


<xe|x_|> 


o> 
g 
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becomes less conservative as the amount of structural damping increases and 
as the amount of nonlinear drag increases. The linear drag factor estimated 
in this manner does not become unconservative (except in the low frequency 
region) for the range of variables studied. 

The amount of conservatism introduced by the use of 

xe 1x |> 

“> 

in linearizing the drag factor, Ky» is seen in Figs. 2.17 through 2.19. In 
these figures, the spectrum for the linear drag solution is seen to be con- 
servative up until the high frequency tail begins. The linearized spectra 
were found to be slightly unconservative, as expected, in the low frequency 
region, but the scale of the figures is such that this unconservatism cannot 


easily be seen. The use of Eq. 2.4, repeated here for convenience, 
“2 


KD : 2 KD Xo max Xg max 
IZ max i (4) (XG ae K ) (x ) (2.4) 
g max 


is also seen to provide a very good estimate of the nonlinear drag spectrum 
in the high frequency region. At the point in the spectrum at which the high 
frequency tail just begins to appear, the intersection of the linear drag 
Spectrum with the estimate of the high frequency tail provided by Eq. 2.4 

is seen to be slightly unconservative. The amount of unconservatism in this 
transition "saddle" to the high frequency tail is reduced as the amount of 
nonlinear drag in the system increases. It is seen from the figures that, 
while linearizing the drag force by the relation, 


“xe |x |> 
op = Ky (a) (2.8) 
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yields conservative results in the mid-frequency region of the spectrum, the 
degree of conservatism is not so great as to be considered too conservative. 
The presence of the nonlinear hydrodynamic drag force has been seen 
to cause a reduction in eee from that computed for the "no-drag" oscillator 
at a given frequency, up until the high frequency tail region is reached. A 
measure of how much of this reduction is accounted for in the computation of 
Te nay for the system with linear drag (linearized by the ground velocity) 
is seen in Fig. 2.20. The ordinate in the figure is the ratio of the reduc- 
tion (from the "no-drag" case) in Z| when the actual nonlinear drag force 
is used to the reduction (again, from the "no-drag" case) in 1a are when the 
drag force has been linearized by the ground velocity terms. In Fig. 2.20, 
values of the ordinate greater than 1] represent unconservative results. The 
unconservative nature of the linearization in the low frequency region and 
the high frequency tail region is noted. While the linearization of low-drag 
systems is seen to account for less than 50 percent of the actual reduction 


in’ | Z| in many cases, the actual reduction for low-drag systems is quite 


ma x 
small itself, as seen in Fig. 2.17. It is further seen from the figure that 
more of the reduction in ae 1S accounted for in the high-drag systems than 
in the low-drag systems. 

Figures 2.17 through 2.20 show a measure of the amount of conserva- 
tism introduced in Ele by using the ground velocity characteristics to 
linearize the drag force when the structural damping ratio, Bas 1S20.02 
For systems with less structural damping, the amount of conservatism will 
be slightly greater, and for systems with more structural damping the amount 


of conservatism will be slightly less. 
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2.4 Linear Hydraulic Damping 


The preceding section of this chapter has demonstrated that the 
drag factor Ky can be linearized in a manner which yields slightly conserva- 


tive values of ve Over a wide range of frequencies in the response 


ax 
spectrum. The linearized drag factor, KOT? can be converted into a linear 
hydraulic damping ratio, By» as shown by Eq. 2.7, which is repeated here for 
convenience. 


- ; ° ee ; 
Zo + 28 wZ + 28 wZ + wZ = “Xo 28, WX 4 (Za) 


Equation 2.7 is seen to have the same form as Eq. 1.12 with the exception of 
the coefficients of the forces on the right side of the equation. Since the 
solution of Eq. 1.12 was shown in Chapter | to be necessary in order to use 
the normal mode method of solution for the response of the multi-degree of 
freedom model, the effects of By on the solution of Eq. 2.7 will now be in- 
vestigated. 

Equation 2.7 is solved numerically for Rae for all four earth- 
quakes with values of B. ranging from 0 to 0.10 and values of By ranging from 
O to 0.10. These ranges of B. and By are considered to be inclusive of all 
reasonable values of these variables for the type of structure considered in 
this thesis. The results of the computations are shown as response spectra 
in Figs. 2.21 through 2.24. These response spectra are seen to have the same 
general shape as spectra for an oscillator without hydraulic damping until 
the high frequency region is reached. It should be noted that a high fre- 
quency tail also appears in the spectra for systems with equivalent linear 


hydraulic damping. The magnitude of this tail can be shown to approach 
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28, 


z MX, a) from above. However, since the equivalent linearization of the 


nonlinear hydrodynamic drag force has already been shown to break down as 
the drag force becomes more dominant in the solution at high frequencies, an 
estimate of (ela in the high frequency tail region of the spectrum for a 
given system can best be made by the use of Eq. 2.4. 

Newmark et al.?2>1%>!* construct an idealized representation of 
the single degree of freedom response spectrum for an oscillator without drag 
forces by using amplication factors for displacement, velocity, and acceler- 
ation. These amplication factors are multiplied by the corresponding maximum 
ground motion quantity in order to obtain an idealized spectrum for an earth- 


» and X If 


quake motion of known (or assumed) values of x 9 max’ 


g max’ xg max 
amplification factors can be developed for the response spectrum of the system 
with equivalent linear hydraulic damping when B. and By are known, and the 

high frequency tail of the spectrum for the actual nonlinear system is esti- 
mated from Eq. 2.4, then an idealized response spectrum can be constructed for 
a given set of conditions for the oscillator with nonlinear hydrodynamic drag 
forces. 

During the course of establishing the amplification factors necessary 
for the construction of idealized response spectra for systems with linear 
hydraulic damping ratios (Eq. 2.7), it was found that, within the range of 
variables studied, the amount of total linear damping, B. + Bp» 1s the impor- 
tant factor controlling the response quantity [alee throughout the spectrum 
except for the high frequency region. Comparisons of ape computed from 


Eq. 2.7 at 31 frequencies from 0.04 cps to 15 cps were made for all four 
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earthquake records for total damping ratios equal to .02, .05, .07, .10, .12, 
and .15 with the damping being distributed between structural and hydraulic 

as shown in Table 2.3. The comparison showed that, up to the beginning of 

the high frequency tail in the spectrum, for a given value of the total damping 
ratio, gee computed for the first distribution of 8. and 8 was within 5 
percent of TA computed for the second distribution of B. and By: Further, 
the great majority of the comparisons revealed a difference of less than 2 
percent between Al ae computed for the two different distributions of damping 


at a given value of the total damping ratio. The ratio 


C= for 8. = A, & = B 


{Z| for 8. = B, B =A 


max h 


was found to be sometimes greater than 1 and sometimes less than 1, indicating 
that at least a small part of the difference may be the result of the conver- 
gence criterion (0.1 percent) used in the numerical integration procedure. 


A comparison of {Z| when 6. = 0.00 and 8, = 0.15 with ree when 


max 
Be = 0.15 and Bh = 0.00 was made for each of the four earthquakes considered. 
It was found that, for frequencies up to the beginning of the high frequency 
tail of the spectrum, the value of ae, computed for the first distribution 
of total damping was not always within 5 percent of the value of sae. com- 
puted for the second distribution of total damping. Most of the comparisons 
differed by less than 10 percent, but a few showed differences on the order 
Of 15 percent. 


The conclusion to be drawn is that, for frequencies up to the be- 


ginning of the high frequency tail region, ee computed from Eq. 2.7 is 
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dependent upon the total linear damping in the system, and not the distribution 
of damping between structural damping and equivalent linear hydraulic damping 
for values of Bh up to and including 0.10. Thus, the use of the amplification 
factors published in Ref. 14, and included in this thesis as Table 2.4, is a 
valid and proper method of constructing a portion of an idealized response 
spectrum for the single degree of freedom oscillator with drag characteristics. 
The scale of the response spectra shown in Figs. 2.21 through 2.24 
is so small that the conclusion is difficult to see from the figures. For 
illustrative purposes, a comparison of response spectra for B, = .07 is shown 
in Fig. 2.25 for the Ferndale earthquake. It is readily seen from this figure 
that, up until the high frequency tail begins to appear, the spectral values 
are a function of the total amount of linear damping in the system, and the 
relative distribution of the total damping between structural and hydraulic 


has little effect. 
2.5 Idealized Elastic Deformation Response Spectrum 


It was mentioned in Section 1.5.1 that the characteristic shape of 
response spectra for earthquake motions has led to the development of amplifi- 
cation factors from which an idealized response spectrum can be constructed 
for a given set of conditions. The idealized spectrum permits an estimate 
to be made of the maximum response of the system. This response quantity 
may then be used as a modal response quantity for a particular mode of a 
multi-degree of freedom system, and the various modal response quantities 
may then be combined in some manner to yield an approximate solution for the 
multi-degree of freedom system as outlined in Chapter 1. 


The principles and concepts developed in Sections 2.2 through 2.4 
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are summarized here for the reader in such a manner as to demonstrate the 
construction of an idealized response spectrum for a single degree of freedom 
system, with nonlinear hydrodynamic drag forces, subjected to a strong-motion 
earthquake of the type used in this thesis. 
To construct an idealized deformation response spectrum for a given 
system: 
1. For the system, compute the drag factor, Ky and the undamped 
circular natural frequency of the system, w. 
2. Establish the structural damping ratio, Bo» to be used in the 
system. For the type of structure discussed in this thesis, 
a structural damping ratio of between 0.005 (when the stress 
is below one-quarter of the yield stress) and 0.05 (when the 
stress is at, or just below, the yield stress) is reasonable 
to assume.’** 


3. Select the values of x for which the 


g max’ XG max” *g max 
spectrum is to be constructed. 


4. Compute the quantity 


“xe | |> 


K 
aD es 





B = 
n 2 Mw 


[It should be noted from Table 1.1 that the magnitude of 


ae x |> 





<Ko- 
g 


1S approximately one-half the magnitude of X a. for the four 


earthquakes shown. 





Se 


Enter Table 2.4 with the value of B. expressed as a percent, 
obtain the displacement amplification factor, multiply the 
displacement amplification factor times Xq 9 and plot the 
result on tripartite logarithmic graph paper as a displacement 
line. This is shown as line 1] on Fig. 2.26. The method of 
linearizing the drag force by using ground velocity character- 
istics has been shown to be unconservative in the low frequency 
region. The use of Be alone means that any hydraulic damping 
effects in the low frequency region are neglected. This is 
reasonable, however, because the small amount of velocity re- 
sponse in the low frequency region results in minor amounts of 
hydraulic damping. Further, as the frequency of the system 
becomes lower and lower, the amount of damping of any type has 
little effect on the maximum deformation. 

Enter Table 2.4 with the amount of total damping (B. + B, as 
computed above) and obtain velocity and acceleration amplifi- 
cation factors. Multiply the velocity amplification factor 
times x and plot the result as a horizontal line. Multiply 


g max 
the acceleration amplification factor times Xo tg and plot the 
result as an acceleration line. These two lines are shown as 
lines 2 and 3 on Fig. 2.26. 


K 
Plot the quantity (2) 2 as a displacement line. This 


xg max 
line is shown as line 4 on Fig. 2.26. 

The envelope formed by the intersection of the four lines plotted 
is the idealized response spectrum for the system and conditions 


given. 
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With the idealized response spectrum constructed as described in 
(1) through (8) above, an estimate of the maximum response of the system is 
obtained by entering the figure with the natural frequency of the system and 
), or pseudo-acceleration (aaliz’ 


reading |Z|___, pseudo-velocity (w|Z| 


max max Fai) 


If the characteristics of the system are changed, a new response spectrum 
Should be constructed. 
Linear interpolation can be used to obtain amplification factors 


for values of damping not shown exactly in the table. 
2.6 Hydrodynamic Parameters 


The hydrodynamic force exerted by the fluid on the mass of the oscil- 
lator as it moves through the fluid is represented by Morison's formula, Eq. 1.2 
as described in Section 1.2. This force expression contains two dimensionless, 
empirical coefficients, Cy and Cy» which are dependent upon the response of 
the oscillator for their magnitude. Many investigators have undertaken 
Studies to determine the value of these coefficients under various conditions. 
The results of some of these investigations are published and discussed in 
Refs. 17 and 18. Morison's investigations” show a variation of the drag co- 
efficient Cy with Reynolds Number, which is related to the relative velocity 
between the mass and the fluid. Keulegan and Carpenter!? have shown a variation 
in both coefficients with a dimensionless period parameter which can be in- 
terpreted as the ratio of the relative displacement between the mass and the 
fluid to the pile diameter for forces against circular piles. Laird?® has 
Shown that the two coefficients may vary greatly for flexible oscillating 


cylinders. Evaluation of the two coefficients for oscillatory motion is further 
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complicated by turbulence associated with eddy formations on either side of 
the mass as it moves back and forth through the fluid.*® This thesis con- 


siders that the drag factor, K,, and the added mass of the oscillator, m , 


D 
remain constant throughout the entire time history of the oscillator motion. 


The values of K, and m_ used in the calculations in this thesis do not pre- 


D 
suppose any specific values for Cy and Cys except as values were assigned 

to these coefficients by Burke* in the formulation of properties of a series 
of offshore towers which are used as sample structures in a later chapter. 
However, the designer of an actual offshore structure, whose physical pro- 
perties and dimensions will eventually be known, must have a means of esti- 
mating the coefficients Cy and Cy with reasonable accuracy in order to analyze 
the response of his structure. Further, investigators who are engaged in 
experimental research on the response of submerged or partially submerged 
bodies to earthquake excitation should have a method of estimating the coef- 
ficients in order that the experimental work may be conducted over a reason- 
able range of variables. Since C, and C, are most commonly considered to be 


J D 


dependent upon the relative motion between the mass and the fluid, X and Xe 


in the case of the mass moving through the fluid at rest, a method is developed 


in this section for estimating the response quantities Ix. | <[x_|>s |x 


max’ Blase 


<|x, |> (hereinafter called hydrodynamic parameters) for strong motion earth- 
quakes of the type considered in this thesis. 


Consider the equations 


MZ + C2 + Ky) 2 + KZ = “MK, - KOL%g (2.6) 


and 
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i + 28.ul + 2Bul + a “iy = 28, (2.7) 


is which the nonlinear drag term Ky has been linearized in the manner dis- 
cussed previously. In the response computation, the quantities Xe and ae 


may be computed at each time point from the relations 


and 


Equation 2.7 is solved numerically for values of 8. and 8, which are reasonable 
to expect for the type of system considered, and again it is found that the 
hydrodynamic parameters 1X. | 


<|x.[>; Ix. | , and <|x.|> are affected by 


max? max 


the total amount of damping in the system and not by the distribution of damp- 
ing between structural and hydraulic. Comparisons were made of the hydro- 
dynamic parameters when the total amount of damping was equal to 2, 5, 7, 

10, 12 and 15 percent of critical with the damping being distributed between 
Be and By as shown in Table 2.3. The difference in a hydrodynamic parameter 
at a given value of total damping ratio is within two to five percent for the 
different distributions of damping between Ee and BA: thus indicating that, 

if the amount of total damping can be estimated, the parameter can be esti- 
mated. The comparison holds true for all four hydrodynamic parameters. 


Spectra for the hydrodynamic parameters Ix. | <|x.|>> [x 


max’ wl 


and <|x.|> are shown in Figs. 2.2/7 through 2.30 for total damping ratios up 
to ten percent. The general shape of the spectra indicates that an idealized 


spectrum may be able to be constructed for a given set of conditions in much 
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the same manner as an idealized deformation response spectrum is constructed. 
Each hydrodynamic parameter spectrum was plotted, a smooth spectrum which 
approximated the actual spectrum was established by the author, and appro- 
priate amplification factors were computed for each of the four earthquakes 
considered. These amplification factors were then averaged in a mean square 
sense to yield the amplification factors listed in Table 2.5. The use of 
these amplification factors to construct idealized spectra for the hydro- 


dynamic parameters x. | <|x.|>; | x » and <|x.|> is explained in the 


max’ Blak 


next section of this chapter. 
2./ Idealized Elastic Hydrodynamic Parameter Response Spectra 


The construction of idealized response spectra for the hydro- 
dynamic parameters discussed in Section 2.6 is very similar to the construc- 


tion of an idealized elastic deformation response spectrum. 


2./.1 Displacement Parameter 


To construct an approximate response spectrum for |x. | (or <|x.|>) 


max 
for a given system: 
1. Determine the quantities Kp Ws Bos and Bh in the same manner 
as described in Section 2.5. 
2. Enter the appropriate section of Table 2.5 with the value of 
B. expressed as a percent and obtain the displacement ampli- 
fication factor. Multiply this amplification factor times Xq aoe 
(or Sea) and plot the result as a displacement line on tri- 


partite logarithmic paper. This line is shown as line 1 on 


Pal e201 lh 





Oa 


To construct an idealized response spectrum for |x 


Ad) 


Enter the appropriate section of Table 2.5 with the amount of 


total damping (8. + 8,) and obtain the velocity amplification 


hn) 
factor. Multiply this amplification factor times X a ele 


<[x,1>) and plot the result as a velocity line as shown by 

line 2 on Fig. 2.31. 

Plot the value EE (or a as a displacement line. This 
is shown as line 3 on Fig. 2.31. 

The envelope formed by the intersection of these three lines 


is an idealized spectrum for |x or <|x.|>). 


olrew ( 


Velocity Parameter 


wiles (or <|X.|>) 


for a given system: 


1. 


Obtain and use displacement and velocity amplification factors 
from the appropriate section of Table 2.5 exactly as described 
above for the construction of an idealized spectrum for a dis- 
placement hydrodynamic parameter. The resulting lines are shown 
as lines 1 and 2 on Fig. 2.32. 

Enter the appropriate section of Table 2.5 with the amount of 
total damping (B. i 8) and obtain the acceleration amplifica- 
tion factor. Multiply this amplification factor times Xo — 
(or SS Pa) and plot the result as an acceleration line. This 
is shown as line 3 on Fig. 2.32. 


max 
by line 4 on Fig. 2.32. 


Plot the value of X (or <[xg|>) as a velocity line as shown 
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4. The envelope formed by the intersection of the four lines is 


an idealized response spectrum for | x (or <|x,|>). 


silniax 
Linear interpolation can be used to obtain amplification factors 


for values of damping not shown exactly in the table. 
Unfortunately, no reasonably consistent ratios were found to exist 
between aieal <|x, I>, leer and other ground motion quantities such as 


Ka x? Ko Bs Rg oe for the four earthquakes considered. Thus, the de- 


Signer of an offshore structure who wishes to estimate the hydrodynamic param- 
eters <|x,|> and <|x.|> will have to obtain the values of = [al <|Xg|>, and 
aK! from processing an actual earthquake record or estimating them using 
the values shown in Table 1.1 as a guide. The range of the ratios between 


IK, |> <[xQ|>s <|KQ/> and x are repeated here for the 


: * 
g max” “g max’ xg max 
convenience of the reader. 


Ratio Lowest Value Highest Value 
<|x_|> 
0.194 O53S5 
g max 
<|x |> 
0.124 0.167 
g max 
<|X |> 
< 0.084 0.185 
g max 


It should be noted that the idealized spectra for hydrodynamic param- 
eters become less descriptive of the oscillator response at very low frequencies, 


about 0.08 cps and less for the four earthquakes used in the computations. This 


<[x.|>> |x 


is to be expected, since the quantities 
E q allen S s |max 


» and <|x.|> 


become quite small for very soft (low frequency) systems. 
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2.8 Inelastic Response Spectra 


In this section, the effects of the nonlinear hydrodynamic drag 
force on the response of the oscillator in the inelastic range are investi- 
gated. Figure 2.33 shows the resistance-deformation relationship of the 
spring in a system with elasto-plastic spring characteristics. The character- 
istics of the elasto-plastic system are an initial elastic region, a constant 
plastic resistance limit, and an elastic unloading region. The unloading 
portion is considered to have the same slope as the initial elastic curve, and 
the unloading is elastic until yielding is reached in the opposite direction. '” 
The yield point deformation is denoted by Ly» and the absolute maximum dis- 


placement of the mass relative to the base is denoted by |Z| The yield 


max- 
resistance is denoted by Ae The ductility factor, u, is defined as 


|Z | 
__ max (2.9) 


i 
y 

and expresses the maximum deformation of the inelastic system in terms of its 
yield point deformation. With up defined in this manner, the maximum deforma- 
tion of the system is Hs the maximum elastic deformation is Zi and the 
maximum inelastic deformation is (u - 1)2,. A variety of calculations of 
the response of inelastic systems has been made and reported in the litera- 
ture.’+~!* The results of the calculations are conveniently displayed in 
the form of response spectra in much the same manner as response spectra 
for elastic systems are displayed. For elasto-plastic systems, response 


spectra can be plotted for constant values of the ductility factor, nu, and 


when the plot is made in terms of the elastic component of the displacement, 
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the accelerations are properly presented.’* In this type of spectrum represen- 
tation, the pseudo-velocity quantity plotted is we Elasto-plastic response 
spectra for systems with 6. = 0.02, KX g max/k = 9-00 (no drag), subjected 

to the El Centro earthquake are shown in Fig. 2.34 for values of u = 1 (elastic 
case) to p = 4. In these spectra, the total deformation of the system can be 


determined by reading the displacement quantity, which is Lys and then using 


the relationship 


Z| max : uy 


In the low frequency range of the spectra shown in Fig. 2.34, the 
total deformation is approximately the same for the elasto-plastic systems 
as for the elastic systems. In the high frequency region, the pseudo-acceler- 
ation is approximately the same for the elasto-plastic systems as for the 
elastic system. In the mid-frequency region, some generalizations have been 
made regarding the amount of reduction in yield level for different values 
of u, but more detailed generalizations are still under study.?? 
Elasto-plastic deformation response spectra for systems with non- 
linear hydrodynamic drag forces subjected to the El Centro earthquake are 
Shown in Figs. 2.35 through 2.37. Again, in the low frequency region, the 
total deformation is approximately the same for the elasto-plastic systems 
as for the elastic system. In the mid-frequency region, the spectra have 
the same general shape as do the elasto-plastic spectra for systems without 
the drag forces, but no quantification of the results has been made. Further 
Study in the direction of establishing rules for the construction of idealized 


elasto-plastic response spectra in this frequency region is clearly indicated. 





48 


In the high frequency region, the elasto-plastic spectra for the systems with 
the lowest value of the drag parameter (0.006) approach approximately the same 
values in the high frequency tail as the elastic spectrum does. However, for 
the systems with the greatest value of the drag parameter (0.50), the high 
frequency tails of the elasto-plastic spectra approach values which are re- 
duced by as much as 50 percent from the elastic case. In other words, there 
is a reduction in the magnitude of the elasto-plastic high frequency tail from 
that obtained in the elastic case, and the amount of the reduction 1s a func- 
tion of the amount of nonlinear hydrodynamic drag in the system. In the sys- 
tems with a small amount of drag, the reduction is very slight, and may not 
merit much consideration, but the reduction for high drag systems can be sig- 
nificant. A study of the numerical results obtained in this investigation 


revealed that, for the spectra with 8. = 0.02, K,x /K = 0.50, values of 
S Dg max 


yu up to 4 and the El Centro earthquake, the factor a=] applied to the 
elastic case high frequency tail yields a satisfactory estimate of the elasto- 
plastic high frequency tail (see Table 2.6). The reader is cautioned that 
this expression for the reduction factor is applicable only to the specific 
systems and earthquake used, and is not to be extended to other conditions. 

In order to compare the effects of the nonlinear hydrodynamic drag 
force on the elasto-plastic response spectra, spectra are plotted for given 
values of structural damping ratio, Bo. and ductility factor, u, for systems 
with varying amounts of drag subjected to the El Centro earthquake. These 
comparisons are shown in Figs. 2.38 through 2.40. These figures show that, 


for a small ductility factor (uy = 1.5), increasing amounts of nonlinear hydro- 


dynamic drag reduce the value of the yield deformation, Zs for a range of 
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frequencies between low frequencies and the beginning of the high frequency 
tail. This means that the maximum response of the system with high drag forces 
is less than the maximum response of the system with lower drag forces, or, 
the drag forces tend to produce a damping effect on the maximum deformation 
of the spring in the mid-frequency region of the spectrum. For example, in 
Fig. 2.38, at a frequency of | cps, and the highest value of the dimensionless 
drag parameter used (0.50), Z is only about 50 percent of the value of Z 
for the system with no drag. However, as the ductility factor u increases, 
the effects of the nonlinear hydrodynamic drag force become less and less 
noticeable in the mid-frequency range. In Fig. 2.40, with u = 4, the value 
of is for the highest drag system is about 90 percent of the value of Z for 
the system with no drag at a frequency of | cps. This is reasonable to expect, 
Since the velocity of the oscillator, Xe can be expected to be smaller during 
plastic deformation of the spring than during elastic deformation of the 
Spring, and the magnitude of the drag force is a function of the oscillator 
velocity. It is also to be noted from Fig. 2.40 that it is possible for Z. 
in the case of the system with drag to exceed slightly the value of Z for 
the system with no drag at the same frequency (note the response values at 
frequencies of 0.8 cps and 3.0 cps) in the mid-frequency range. This can re- 
Sult since the maximum response of the elasto-plastic system depends on the 
entire history of the deformation and the energy absorbed in the spring during 
plastic deformation as well as the energy dissipated by the damping effects 
of the fluid surrounding the oscillator. 

The significant features of the elasto-plastic response spectra 


for systems with nonlinear hydrodynamic drag are: 
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The drag does not have much effect on the spectrum in the low 
frequency region; 

A high frequency tail does exist, and large amounts of drag 

can cause significant reductions in a for increasing values 

of the ductility factor (less significant reductions in Z. with 
increasing values of yw for low drag systems); 

The damping effects of the drag in the mid-frequency region 

of the spectrum become less pronounced as the ductility factor 


increases, even for the systems with high drag. 
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CHAPTER 3 
MULTI-DEGREE OF FREEDOM SYSTEMS 
3.1 General 


In this chapter, the principles developed and discussed in Chapter 
2 for the single degree of freedom oscillator are extended to the multi- 
degree of freedom system. The nonlinear drag forces are linearized in such 
a manner as to insure the existence of classical normal modes in the equa- 
tions of motion in order that normal mode superposition may be used as an 
approximate solution to the problem. The method of linearization is shown 
to be accurate for a series of five degree of freedom systems which are 
representative of the type of structure being investigated. A means of es- 
timating the amount of equivalent linear hydraulic damping, using ground 
motion velocity characteristics, is developed and is shown to yield good re- 
Sults. The construction of an idealized response spectrum for each mode of 
the multi-degree of freedom system is discussed and described, and sample 


problems are solved by the superposition of normal mode responses. 
3.2 Approximate Solution of the Nonlinear Equations of Motion 


In this section, an approximate method of solving the nonlinear 
equations of motion for the multi-degree of freedom system is developed and 
Shown to yield satisfactory results. Various features of the solution, such 
as the formulation of the structural damping matrix, the linearization of 
the drag forces, the superposition of modal responses, and the response of 


very high frequency modes, are discussed. The transformation of the nonlinear 
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equations of motion into uncoupled linear modal equations, which was first 
shown in Chapter 1, is repeated in this section for the convenience of the 


reader. 


3.2.1 Equations of Motion 


The equations of motion for the lumped parameter model shown in 
Fig. 1.1(b) are written in matrix form as Eq. 1.3 and are repeated here to 


form the basis of the analytical discussion in this section. 
PM JCZ) + [C.1€2) + EK]{Z} = -P MIG} - DK fg 1k |} (3-1) 


No structural dampers are shown on the model, but the damping matrix LC] 
can take a variety of forms as will be discussed in later paragraphs. For 
a given structural system and known earthquake record, response quantities 
of interest may be obtained by the numerical integration of Eq. 3.1, but 
this procedure is costly and requires the use of a digital computer. If 
Eq. 3.1] can be satisfactorily approximated by a system of simultaneous 
linear differential equations of motion which possess classical normal 
modes, then a solution can be obtained by the superposition of normal mode 
solutions. A furtner approximation can be made by using response spectra 
to obtain the maximum response in each mode and then combining the modal 
maxima in some fashion, very much like the approximate modal techniques, 
using response spectra, which are reported in the literature!?~!% for land- 
based structures. 

In the normal mode method, the vector of a response quantity ex- 


pressed in one coordinate system, such as the relative displacement between 
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each node and the base, {Z}, is transformed into the "normal" coordinate sys-~ 


tem by the relation 
{Z} = [o]{yY} (362) 


If the nonlinear drag force were not present, Eq. 3.2 and the ortho- 
gonality relationship between [o] and [-M.], [K], and [C.] could be used to 
transform Eq. 3.1 into nm (number of modes) uncoupled single degree of free- 
dom equations in the normal coordinate system. If the nonlinear drag force 
can be converted into an equivalent linear drag force, then the resulting 
linear differential equations of motion can be transformed into uncoupled 
single degree of freedom equations of motion in the normal coordinate system 
as follows: Consider that [~Kp~J (x. |x, |} can be transformed into an equiva- 


lent linear force, [Ky ~ {x.] j.e., 
DK ytx[x.[} + OK) J {x} (8.5), 
Substitution of Eq. 3.3 into Eq. 3.1 yields, 
PD MJ{Z} + (C12) + [KZ = -OMAGK} - Pky Tk} (3.4) 
Since {2} = od - {xg}; Eq. 3.4 can be rewritten as 
DMZ} + [C.1{Z} + Pky HZ} + [KZ 
= -D MEK} - PK tx,} (ae) 


Substitute Eq. 3.2 into Eq. 3.5 and premultiply the result by co]! 
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Co} OM Ico) + Col" Ec Ilo] t¥} + Col Ky Ilo]t¥} + Col" EKILOI(¥} 
= -[o] DM If%)} - Lo] Pky I) (3.6) 


Recall the following definitions from Chapter I: 


CMW] = Co] EMCO] = a diagonal modal mass matrix 

ig - Co] '[¢_IL¢] = a diagonal modal structural damping 
matrix. The form of [C.] which permits 
this is discussed later. 

PK = [ro] CKICo] = a diagonal modal stiffness matrix 


a | 
[Ky J = (ol Kp) JLo] a modal linear hydraulic damping 
matrix which 7s not necessarily 


diagonal 


* = * 
Bae [~M ~] The ~] a matrix of circular natural fre- 


quencies of the system 
Substitution of these definitions into Eq. 3.6 yields 
* oe * ° Ny ° * 
ee ty) +t [~C ~]{¥} 1s LK, J{Y} ele eee 
> fol Pmt) - bok, 33 (3.7) 
g DL g 


At this point, the formulation of the structural damping matrix, [C.J], and 


the linearization of the drag factor, [~Kp~J, will be discussed. 


3.2.2 Structural Damping 


If the structural damping matrix [C.] is proportional to the mass 





25) 


* 
matrix, then [C~] will be diagonal since the natural mode shapes are ortho- 
gonal with the mass matrix. This formulation of [C.] results in modal damping 


ratios which decrease with increasing natural frequency. This is shown as 





fol lows: 
Let, 
[c] = etmy 
then, 
T T * * 
Col'(c Jo] = éfel'CMJte] = et = Ce (3.8) 
Define 
* =] * 
CM oc = C2ewJ = et] (3.9) 
then, 
eee 
B = (35g) 
S. 20), 
where 
B. = structural damping ratio in the ith mode 
F 
W.. = circular natural frequency of the ith mode 


In a similar manner, it can be shown that, if the structural damping 
matrix [C_] is proportional to the stiffness matrix [K], then [CW] is a 
diagonal matrix in which the modal damping ratios Bs. increase with increasing 
frequency. 

The elements of the structural damping matrix [CJ], can be selected 
to yield the same structural damping ratio in each mode. It is common prac- 


tice to assume the same structural damping ratio in all modes in the dynamic 
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analysis of damped multi-degree of freedom systems.'* If the modal matrix 


[>] is normalized such that 


fol EMJEs] = EIQ — (a unit matrix) (3.11) 
Then 

[CJ = CMJlelCc.Itel' tm (3.12) 
And, since 

Pe = Deb 280 (3.13) 


Equation 3.12 can be written as 


[oJ = CMIlelCM Jo 268.0 Jel tM (3.14) 


Thus, it is seen that [C.] can be formulated in such a manner as 
to yield modal structural damping ratios which are constant, which increase 
with increasing mode number, or which decrease with increasing mode number. 
In the computations reported in this thesis, the structural damping matrix 


[C.] is computed in such a way that B. is the same in all modes. 


3.2.3 Linearization of the Drag Force 


2? and Malhotra and Penzien® linearize the nonlinear drag 


Foster, 
forces, [Kp Hix, |x, |} by replacing them with a velocity dependent linear 
force, CK ee, at each node which performs the same amount of work during 
a period of time as does the nonlinear force at that node. Foster then 


recognizes that the resulting [~kK,, ~J matrix is not orthogonal with the mode 


Shapes and proceeds to solve the equations of motion by a nonclassical mode 
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approach using complex eigenvalues and eigenvectors. Malhotra and Penzien 
transform (“Kp I into [Ko] and then force [Ky] to be diagonal by means 
of an iterative process that replaces the nondiagonal matrix with a diagonal 
matrix which dissipates the same amount of energy over a period of time. 

The hydraulic drag forces, whether linear or nonlinear, may be con- 
sidered to act as dampers, that is, they cause energy to be dissipated as 
the system oscillates through the fluid. In the approximate solution being 
developed in this section, the nonlinear dampers in the multi-degree of 
freedom system are replaced with linear dampers (except for very high fre- 
quency modes which are discussed in a later section) which dissipate the 
same amount of energy over the total system, but do not necessarily preserve 
the dissipation of energy at each node. Further, the resulting “equivalent” 
linear damper at each node is made to be proportional to the lumped mass of 
that node in order that the modal equations of motion will be uncoupled. For 
an offshore structure with a portion of the structure above the water line 
as shown in Fig. 1.1(a), this procedure results in a linear hydraulic damper 
being provided at the level of the top mass where, in fact, no nonlinear 
hydraulic damper really exists. With this exception, the proportioning of 
the linear hydraulic dampers in the same ratio as the mass is distributed is 
a reasonable approach for the type of structure being discussed, since the 
nonlinear damper is proportional to the projected area of the structure lumped 
at the node, and the magnitude of the mass at a given node will also be gen- 
erally proportional to the projected area. The method of linearizing is as 


follows: 





a nn t 
Energy porate ae | (kK, SI eee 
nonlinear dampers 71 } 53 S. S S 
y : 
= t Rk SS ise fe 
ie au) 3 
Energy dissipated by nn ¢t ; , 
equivalent linear = } | (Ky x. )(x, decit 
dampers te Id, 3 J 
nn : 
= t ) K Nee 
in 5 
where 
nn = number of nodes in the system 
j= jth node in the system 
Equating the energy dissipated, 
nn nn 
Ko. elke 1 = DK, <P 
Je jJ WA Vel jJ J 
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seo) 


(3G) 


oll) 


It is noted that only the total energy dissipated is preserved and the dis- 


tribution of the energy dissipation among nodes is not preserved. Now, by 


forcing the linear dampers to be proportional to the mass matrix, 
[Ky~] = ofM. 


and substituting Eq. 3.18 into Eq. 3.17 


(3.18) 





ae) 


nn 2 
deh ea J 
— anu aie a (2en9) 
Py 3 
ii 
a =~ 


AY) 
Recall that [Ky J was defined as Col’ Ky ~IL91. Substicutingsla.. 3.18 into 
this definition, 
AV) T * 
[KJ = ofl’ CMJ0e] = afm (3.20) 
Substituting Eq. 3.20 into Eq. 3.7 yields 


cM ]{¥} + rc.Jt¥} + ofM J(Y} + BK 0} 


= fol EM II-%, - ak) (3.21) 


Equation 3.2] represents nm single degree of freedom uncoupled equations of 
motion in the normal coordinate system. 

Before proceeding, the method of linearizing described by Eqs. 3.18 
and 3.19 will be tested for accuracy in computing the response of multi-degree 
of freedom systems to earthquake excitation. A series of five degree of free- 
dom systems, representative of offshore towers in mass, stiffness, and drag 
distribution, is developed in order to test the linearization technique. 

Seven structural systems, consisting of different combinations of 
mass and stiffness distribution, are considered. Each system is further 
subdivided into one or more distributions of the drag factor, Kp» at each node. 
The characteristics of these systems are shown in Table 3.1. In the table, 
the values shown for the mass at each node are considered to be the total 


mass at that node, and include the added mass term. The different 
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distributions of drag factor are labeled L, M, and H in the table, and the 
column under each of these headings represents the relative magnitude of Ky 
at each node for the particular distribution shown. For example, drag factor 
distribution M represents the distribution of Kp at each node as 0, 1, 1, 2; 
3 at nodes | through 5 respectively. The values of »KA/ 2m Shown in the 

table for the distributions of drag factor are of sufficient range to span 
reasonable values of this quantity for actual offshore towers of the type 
being investigated. The quantity UK / zm for the four towers reported by 
Burke and Tighe” range from a low of 0.0034 inches~! to a high of 0.0054 
inches”! The relative magnitudes of the undamped natural frequencies of 

the system, normalized to the fundamental frequency, are as shown. In the 
computations which are discussed in the following paragraphs, one of the 
natural frequencies is set to a desired value and the fundamental frequency 
is computed from the ratio of natural frequencies. The value of kK, 1S com- 
puted from f, and the quantity f,/vk,/m, » then the remaining spring stiffnesses 


are computed from k, and the distribution of k. through k The stiffness 


1 2 oe 
matrix [K] for the entire system is then assembled and is a tri-diagonal 
matrix for the systems shown. 

Each of the systems shown in Table 3.1 is subjected to the El Centro 
earthquake and response quantities are computed by the numerical integration 
of the equations of motion using the criteria for convergence and time inter- 
val described in Chapter 1. For each system, computations of maximum re- 
Sponse are made for two values of B. in each mode, 0.01 and 0.05, and three 


frequency ranges. The low frequency range, range number 1, is established 


by setting the fifth natural frequency equal to the frequency at which Xq max 
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and x intersect when plotted on tripartite logarithmic paper. The 


g max 
mid-frequency range, range number 2, is established by setting the funda- 
mental frequency equal to the frequency at which Xa a and X Ls inter- 
sect. The high frequency range, range number 3, is established by setting 
the fundamental frequency equal to the frequency at which X eee and Xa ee, 
intersect when plotted on tripartite logarithmic paper. The maximum dis- 
placement of each node relative to the base, UZ ees is computed. The 
results of the computations are shown in Tables 3.2 through 3.27. The 

column in the tables titled "No Drag" shows the results of the computations 
when [~K,~J is set to zero. The column in the tables titled "Actual Drag" 
Shows the results of the computations when [ Ky~J 1s computed from the values 
shown in Table 3.1, thus nonlinear hydrodynamic drag forces, T Kp Hix, 1X, |} 
are acting on the structure. In the computations for the "Actual Drag" 
systems, the time averages (<x |x |>) and {<xo>} are collected. These time 
averages are then used to linearize the drag factor and set it to be pro- 
portional to the mass matrix as prescribed by Eqs. 3.18 and 3.19. Equation 
3.4, using the linear drag factor matrix just computed, is then integrated 


numerically to obtain the response quantity {|Z| The results of this 


max? 
integration are labeled "Case 1" in the tables of results. The column 
labeled "Case 1/Actual Drag" in Tables 3.2 through 3.27 represents the com- 
parison of the solution of Eq. 3.1 with the solution of Eq. 3.4. A study 
of these results reveals that the linearization technique described by Eqs. 
3.18 and 3.19 is extremely accurate. The difference between the linear and 
the nonlinear results is less than 6 percent (with one exception--system 1M, 
B. = 0.01, node 3, mid-frequency range) and is less than 2 percent in the 


majority of the cases. 
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The maximum distortion between two adjacent nodes is also computed 
for both the linear and the nonlinear equations of motion. The results of 
the computations for the linear case in which the linearization is carried 
out according to Eqs 3.18 and 3.19 are not shown in this thesis, but the 
linear solution is as close to the nonlinear solution as for the response 
quantity 2 ee shown in the tables. Thus, it is concluded that the meth- 
od of linearizing the drag force described by Eqs. 3.18 and 3.19 yields 
satisfactory results. 

The primary difficulty in using Eqs. 3.18 and 3.19 to linearize the 
drag force lies in the ability to determine the quantities (<xe |x, |>} and 
(<x6>} a priori. If these quantities can be estimated in some fashion, and 
if the estimate can be shown to be accurate in certain frequency ranges, 
then Eq. 3.1 can be linearized into Eq. 3.4, and Eq. 3.4 can be transformed 
into uncoupled, modal, single-degree of freedom equations in the normal co- 
ordinate system which can be solved by spectrum techniques. It was demon- 
Strated in Chapter 2 that, for the single degree of freedom oscillator, the 
quantities cK IX > and <> are conservative estimates of the oscillator 
velocity in all but the low frequency region of the spectrum. Since the 
major effects of hydraulic damping occur in the velocity amplification re- 
gion of the spectrum, it seems reasonable to estimate the quantities 
(<ke|x_[>} and {<xo>} as (<ko 1X9 [>) and (<Ke>} respectively. Using these 


quantities, Eq. 3.19 becomes 


nn om 
i we, <Xq/Xql? 
af Jel JJ 
a = ae Se (3.22) 
) exe > 
= jj g 


C3, 
—| 
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A series of computations, labeled "Case 2" in Tables 3.2 through 3.27, 1s 
made in which the nonlinear drag factor is linearized in the manner pre- 
scribed by Eqs. 3.18 and 3.22. A study of the column labeled "Case 2/Actual 
Drag" shows that the linearization described by Eqs. 3.18 and 3.22 is, in 
general, conservative in the mid and high frequency ranges and is generally 
unconservative in the low frequency range. The degree of conservatism in 
the mid-frequency range is seen to be less than 5 percent in most cases, 
although one case is conservative by 12 percent. In the high frequency 
range, the results are seen to be conservative by as much as 20 percent in 
some cases, especially in the high drag systems, but many of the results are 
conservative by less than 10 percent in the high frequency range. The degree 
of unconservatism in the low frequency range is seen to be less than 10 per- 
cent in most cases, but can range as high as 23 percent. This unconservatism 
in the low frequency range is to be expected, since very soft (low frequency) 
systems can be expected to exhibit little velocity response at each node, 
and so the use of the ground velocity characteristics overestimates the 
velocity of each node. One further observation to be made from a study of 
the numerical results is that the use of the ground velocity to linearize 
the drag force is generally more conservative as the response is reduced 
to a greater degree by the actual nonlinear drag force. In those cases in 
which the Case 2 linearization is unconservative in frequency ranges 2 and 
3, the degree of unconservatism is generally very slight (less than 1 percent). 
The maximum deformation between adjacent nodes is also computed 
and a comparison of the linear drag response (linearized by the ground velocity 


characteristics) with the nonlinear drag response is shown in Figs. 3.1 
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through 3.13. These figures show that the maximum spring deformation com- 
puted when the drag force is linearized by the ground velocity characteristics 
is also conservative in the mid and high frequency ranges and generally un- 
conservative in the low frequency range. It is concluded from these compu- 
tations that the use of Eqs 3.18 and 3.22 is a reasonable way to linearize 

the drag force and yields slightly conservative results except in the low 
frequency range. The conservatism introduced by the linearization is seen 
from the figures to be generally less than 10 percent in the mid and high 
frequency ranges, with some of the results in the high frequency range being 
conservative by as much as 20 percent. The results of the low frequency range 
calculations are generally unconservative by less than 20 percent, but can be 
unconservative by as much as 25 percent. 

The effect of structural damping on the response is quite apparent 
from the tables and figures. For example, the maximum relative displacement 
between the bottom node (mass number 5) and the ground for system 4M in fre- 
quency range 2 is reduced approximately 17 percent by the hydraulic drag when 
the structural damping ratio is 0.0] in each mode (see Table 3.14). Yet the 
response 1s reduced by less than 2 percent by the hydraulic drag when the 
Structural damping ratio is 0.05 in each mode (see Table 3.15). Thus, the 
selection of structural damping ratios to use in the response computations 


is seen to be very important. 


3.2.4 Modal Equations in Normal Coordinates 


The use of Eqs. 3.18 and 3.22 to linearize the drag force has been 


Shown to yield generally conservative results in the mid and high frequency 
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ranges, and unconservative results in the low frequency range. Since the 
linearization prescribed by Eqs. 3.18 and 3.22 insures that the equations of 
motion for the linearized system (Eq. 3.5) will possess classical normal modes, 
the discussion of the superposition of normal modes is resumed. Recall] the 


uncoupled equations of motion in the normal coordinate system. 
* ee * : * ° * 
eee ay} + [~C.~}{Y} PeceMesluyier [eK ay} 
; [o]' DM II-X, S aX} (3.21) 
Premultiplying by [M~J]~! yields 
(¥} + P26 wt} + 028,01} + Doty} 
_ ||) - : 
= [M-~] [6] CMAI{-x, - OX} (3323) 
where 
af-I.] = [28 w] (3324) 
The terms on the right side of Eq. 3.23 can be written as 
e o-lirj4k < ° _ Fe ° 
CM~] [9] DMJ{-X, - OX a : {PF}(-X, - OX) (3.25) 


where 


{PF} = MJT re) hugo (3.26) 


and PF is the participation factor for the ith mode. Substituting Eq. 3.26 
mto EQ. 3.23 


(¥} + P28 wt} + 0260 I{¥} + Duty} 


= {PF} (-X, - OX) Ko a27) 
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For the ith mode, then, 

V ; : 2 7 fs 

Y, + Soa iP oi aaa + wei, = PF AX, OX ) (3.28) 
Since the left side of Eq. 3.28 is linear with respect to the unknown Y, the 
maximum modal response in the ith mode can be obtained by determining the 


maximum response of the single degree of freedom equation, 


y A 2B, 05 iY; < 28 a 4 we Ys = ae - OX (3.29) 
and multiplying this response times the participation factor for the ith mode, 
PF... The nonlinear equations of motion for the multi-degree of freedom system, 
Eq. 3.1, have now been transformed into a set of linear simultaneous differ- 
ential equations (Eq. 3.5) by the use of Eqs. 3.18 and 3.22. The linear 
Simultaneous differential equations are transformed into nm uncoupled, single 
degree of freedom modal equations of the form represented by Eq. 3.28 in the 
normal coordinate system Y by the use of Eq. 3.2. The total response of the 


system is the result of a combination of the various modal responses represent- 


Bemby Eq. 3.28. 


3.2.5 High Frequency Modes 


The existence of a high frequency tail in the response spectrum 
for a single degree of freedom oscillator was demonstrated and discussed in 
Chapter 2. Further, the linearization of the nonlinear hydrodynamic drag 
force was shown to break down when this high frequency tail began to appear. 
A method for determining the parameters necessary to estimate the magnitude 


of this tail for each mode in the multi-degree of freedom system is now 
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developed. Consider the equations of motion for an undamped multi-degree of 


freedom system 
CMJ{Z} + (KHZ = -D MEK} - Dk, | x, 13 (3.30) 
If the system is so stiff that {x} x {xg}, Eq. 3.30 can be written as 
MJ{Z} + [K]{Z} = -DMJ{%} - DK {xg1x4|} (eks1) 


Transform Eq. 3.31 into the normal coordinate system by the use of Eq. 3.2 


and premultiply by ri 
Coy EM Jto1¥} + Col EKILo1¢¥3 
= -[o] TDM K, - Le] ky JOH KI 1) (3.32) 


Or 
* oe * 
[-M ~J{¥} + [7K ~]{Y} 
= -Lo]'DMICK, - Co] OK) 1X,1) (3.33) 
Equation 3.33 represents nm single degree of freedom uncoupled modal equations 


in normal coordinates for the high frequency multi-degree of freedom system. 


Premultiply Eq. 3.33 by cm 7, then 


x |) (3.34) 


{¥} + PutJ}t¥} = -{PFIK - (K-33 
W) = - Xa - {ky (x, ; 





where 


{ky} = DMD ' tel DK) (3.35) 


For the ith mode then, 
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aN 


* 
¥. + w°.Y, = PF x - otk ky) (3.36) 
1 We 4 1 g PF. GaaG ; 





Again, since the left side of Eq. 3.36 is linear with respect to the unknown 
Y, the maximum modal response in the ith mode can be obtained by determining 


the maximum response of 


aN 


* 
D. 
oe 2 = =e] i 1 e 

Y, FOG, = xg F. xg 





| (3.37) 


9 


and multiplying this response times the participation factor for the ith mode. 


It was shown in Chapter 2 that the quantity LYeuley can be estimated for 
eq. 3.3/ as 
K 
A 
PF. 
- noe 2 
anaes: | Soames) (Xora) (3.38) 
meen. 
aa 
and 
ae - AGB? s be) (3.39) 


3.2.6 Superposition of Modal Responses 


The maximum response in any given normal mode of the multi-degree 
of freedom system can be computed from Eq. 3.28. However, the maximum re- 
Sponse in one mode does not necessarily occur at the same time as the maximum 
response in another mode, and the time of occurrence of the maximum response 
in the time history of the modal response is not normally known. The actual 
combination of modal responses in a multi-degree of freedom system is very 


complex and few generalities can be made about the combination. 
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An approximation to the actual maximum response of a multi-degree 
of freedom structure can be made by determining the maximum response in each 
mode, and then combining the maximum modal responses in some fashion. An 
upper bound on the maximum response of the structure can be established by 
assuming that the maximum modal responses occur at the same time in each 
mode and are additive in nature. In other words, the summation of the ab- 
solute values of the maximum response in each mode yields an upper bound to 
the response of the structure. Another commonly used method of combining the 
maximum modal responses is to take the square root of the sum of the squares 
of the modal responses. Other methods of combining modal responses are some- 
times used.** The two most commonly used methods are the summation of the 
absolute values of the modal maxima (ZABS) and the square root of the sum 
of the squares of the modal maxima (SRSS). These two methods are the only 
Ones used in this thesis for illustrative purposes. Tables 3.28 through 
3.33 show the results of the superposition of maximum modal responses applied 
to three of the five degree of freedom systems. The drag force is linearized 
in accordance with Eqs. 3.18 and 3.22 (ground velocity linearization), and 
the amount of linear hydraulic drag BL In each mode is computed from Eq. 3.24. 
The maximum modal response in each mode is computed from the response of the 
Single degree of freedom modal equation (Eq. 3.29), and multiplied times the 
participation factor for the mode. The LABS combination is seen to be higher 
than the true response (as should be the case) and the SRSS combination is 
generally less than the true response. No generalizations as to the relative 
accuracy for different distributions of mass, stiffness, and drag factor can 


be made, however.’* Thus, it is seen that, if the maximum responses of the 
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single degree of freedom modal equations, Eq. 3.28, can be computed, the maxi- 
mum response of the structure can be approximated by a combination of the 


modal maxima. 
3.3 Idealized Elastic Deformation Response Spectra 


It was seen in the preceding section that the maximum response of 
a multi-degree of freedom system can be approximated by combining maximum 
modal responses in some fashion. If the maximum modal responses can be esti- 
mated in a way that accounts for the effects of the nonlinear hydrodynamic 
drag forces in the system, a rapid approximation of the response of the sys- 
tem can be obtained. The principles developed in Chapter 2 for the single 
degree of freedom oscillator and the method of linearizing the drag forces 
developed in this chapter can be combined to construct an idealized elastic 
deformation response spectrum for each mode of the multi-degree of freedom 
System subjected to strong motion earthquake excitation. The procedure for 
the construction of the modal elastic deformation response spectra is now set 
forth. 
1. Compute the natural frequencies, [~w ], and mode shapes, [9], 
for the system. 
2. Establish the structural damping ratio, Bs.» for each mode. 
Va Wiese of Bs. between 0.005 and 0.05 are reasonable to assume 


for the type of structure being investigated in this thesis. 





3. Determine x » xX me » and 
g max g max g max 
°? & 
<x 1x | > 
°? 


<a> 
g 
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for the earthquake to be used. If the spectra are to be con- 
structed for a hypothetical or "design" earthquake, one-half 
of the value of X aie is a reasonable value to assume for the 
quantity 


“xe |x |> 
oes ae 
eT Pe 
*g 
* * 
Compute M.., Ky, PF., a, Bh.» and the high frequency tail 
11 j 1 j 


estimate from the following expressions: 


cM] = Col fm Jf0] 





{Ky} = CM J! tel Dk, I (3.35) 
2 * J-1-,417 
{PF} = [M-~J [o] [DMI{I} (31826) 
nn 
) KD. <xe x |> 
ee ee a3 (3.22) 
y M Bait 
q=1 JJ 
[280] = af I (3.24) 
* 
D, 
High Frequency Tail PF. r \2 
Estimate for the ith = | ——-;— X (3.38) 
Mode we i. Sige 


With the quantities w. 


»» Bo.» Bh.» and the high frequency tail 
11 7] j 


estimate, construct the idealized elastic deformation response 
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Spectrum for the ith mode exactly as described in Section 2.5 
for the single degree of freedom system. 

It is noted that the idealized response spectrum constructed as 
outlined above accounts for the fluid-structure interaction in a slightly 
conservative manner. When an idealized spectrum has been drawn for each mode 
of the system, the maximum modal responses are obtained from the spectra, 
multiplied times the appropriate modal participation factor, and combined in 


some manner to yield an approximate maximum response of the system. 
3.4 Hydrodynamic Parameters 


An approximate method for computing hydrodynamic parameters for a 
multi-degree of freedom system by the superposition of normal mode responses 
is developed in this section. As explained in Chapter 2, these hydrodynamic 
parameters are absolute displacement and velocity responses of the mass (each 
node in the multi-degree of freedom system), whose magnitudes are desirable 
to know in order that the empirical drag and inertia coefficients, Cy and Cy, 
may be estimated in a reasonable way. The hydrodynamic parameter spectrum 
amplification factors developed in Chapter 2 may be used to construct idealized 
Spectra from which the response quantity of interest in each mode can be esti- 
mated. 


Consider again the linear equations of motion, 
DM J{Z} + [C.1{Z} + [K]{Z} = -DMIGK} - Pky The} (3.4) 


Making the substitution {Z} = {x} - {xg}, Eq. 3.4 can be written as 
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PM J{k.} + [C ]{x.} + DK) ~1{x.} + [kK] {x} 
= [K]{x,} + [C.]tx,} (3.40) 


If the velocity force, [C.1{x,}; can be dropped from Eq. 3.40 for values 

of [c.] up to some magnitude without unreasonable loss of accuracy in the 
solution (this is discussed in subsequent paragraphs of the section), and 
the drag factor matrix [kK ~I is made to be proportional to the mass matrix 


[-M.], then Eq. 3.40 can be written as 
PD MACKF + [C ]ix.} + of MI{x.} + [kK] {x} = [K]{x,} (3.41) 


In order to transform Eq. 3.41 into a "normal" coordinate system, the trans- 


formation 
{xt = Lo]{Qh (3.42) 


is made. Substituting 3.42 into Eq. 3.41 and premultiplying the result by 
i 
Lo]. 


Pw JQ} + DC IQ) + of JQ) + DK tO) 
= Lo] Ek] x4) (3.43) 


Define 


fy} = Col EKICI} (3.44) 
and Eq. 3.43 can be written as 


CM HO + CCSD + ol tO + OKI = Ox, 
ioe45) 
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Equation 3.45 represents nm uncoupled single degree of freedom modal equations 


of motion. The equation of motion for the ith mode is 


* fe * ° fh * * _ AG 
ere sp ti i Va%y sili 
or, 
d, + 28, wid, + 28, 0,50, +020, = (ete yu%x.) (3.47) 
i es ha ie ieee 04 mere ii”g 
liza 


Since the left side of Eq. 3.4/7 is linear with respect to the unknown Q, the 
response quantities of interest in the ith mode can be obtained by computing 
the desired response quantities from 


dat 28 Ww i 28, wq or Ae = WX (3.48) 


and multiplying the result times the ith mode participation factor for hydro- 
dynamic parameters. The participation factor for hydrodynamic parameters 


can be seen from Eas. 3.44 and 3.4/7 to be 
D * s-]- 4 
{HDPF} = [“K ~] [o] [K]{T} (3.49) 


Once the modal quantities are obtained in the Q coordinate system, Eq. 3.42 
can be used to obtain an approximation of the actual hydrodynamic parameter 
of interest at each node in much the same manner as an estimate of the maxi- 
mum deformation of each node is obtained by modal superposition. 

To answer the first question raised in this section, the suitability 
of dropping the velocity force, [C.]{x,}, from Eq. 3.40, response calculations 


were made for the following two equations of motion: 
DM {x} + [C.]{x J * (“Kn I{x.3 + LK] {x} 


= [K] {xg} + [C.1{x9} (3.40) 
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and 


PM J{K.} + [C]tx} + DK tx.) + CK] {x,} = [K]}{x} (3.50) 


The quantities {|x bs {<|x.|>}, {[x.| 5,22 and {<|x,|>} for systems 1, 


s|max 
4, and 7 in the mid-frequency range were calculated for the El] Centro earth- 
quake from both Eqs. 3.40 and 3.50. The damping matrices [C.] used in the 
series of calculations were such that they yielded damping ratios of .02, 
.05, and .10 in each mode and [~Ky) J was set to zero. For each structural 
system and value of [c.], the hydrodynamic parameters computed from Eq. 3.40 
were within 5 percent of the corresponding response quantities computed from 
Eq. 3.50. The results of the calculations for system 4 and a damping matrix 
which yields .02 in each mode are shown in Table 3.34. In another set of 
calculations, the damping matrices [C.J were made to be proportional to the 
mass matrix and the numerical values of [C.J were such that the first mode 
damping ratios were .02, .05, and .10. Again, for each structural system 

and value of [c.]. the response quantities computed from Eq. 3.40 were within 
5 percent of the corresponding response quantities computed from Eq. 3.50. 
The numerical results were such that no clear trend resulted that indicated 
which equation yielded consistently higher or lower results. The hydrodynamic 
Parameters were eonsumea for system 4M in the mid-frequency range, with a 
structural damping matrix which yielded modal damping ratios of .02 in each 
mode, for the El Centro earthquake. The drag factor matrix was linearized by 
using the ground velocity characteristics and made to be proportional to the 
mass matrix. In one set of computations, the velocity force [C.1{x,} in Eq. 


3.40 was included and in another set of computations, the velocity force 
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[C.}{x,} was neglected. The two sets of numerical results differed by less 
than 5 percent. 

From these calculations, it is concluded that the velocity force 
[C.]{x,)} can be dropped from Eq. 3.40 without sacrificing the accuracy of 
the results by more than about 5 percent, for values of [C.J which yield 
modal damping ratios of up to about .10. 

In order to show the similarity of form between Eq. 3.48 and the 
linear equation of motion for the single degree of freedom oscillator dis- 
cussed in Chapter 2, it is instructive to recall Eq. 2.7, repeated here for 


convenience. 
ee Ld i 2 _ ve e 
Hate 2B wl + 28) wZ +wZ = “Ko - ote (351 ) 


Since Z = At Xa Eq. 3.51 can be written 


a : ° Z ee ° 

x. + 2B UX. + 2B, WX. twxX, = w Xo + 28, WX (3.52) 
It was seen in Chapter 2 that, for values of Bs otal up to about .10, the 
hydrodynamic parameters computed from Eq. 3.51 were dependent upon the amount 


of total damping in the system (6. + 8) and not on the distribution of the 


hn? 
damping between structural and hydraulic. If all of the damping in Eq. 3.5] 
1s considered to be structural and all of the damping in Eq. 3.52 is considered 
to be hydraulic, the velocity force on the right side of both equations can 


be dropped, and Eqs. 3.51 and 3.52 become 


7 4 28 wi ae -K, (3.51a) 


and 


- ° 2 2 
et 28, WX. See OUX.. “Me Xe (3.52a) 
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In other words, for values of B, up to about .10, the hydrodynamic parameters 


| x <|x.|>> | x , and <|x. [> computed from Eq. 3.5la are within 2 


Blmax’ Mita 


to 4 percent of those which would be computed from Eq. 3.52a for the same 
system and base disturbance. 

It is seen that Eq. 3.48 is of the same form as Eq. 3.52a, so if 
the hydrodynamic parameters can be calculated, or estimated (as from an ideal- 
ized response spectrum) for Eq. 3.52a with satisfactory accuracy, then modal 


hydrodynamic parameters Id ax? <|q|>, |q|_.., and <]q|> can also be calculated 


max 


or estimated. Once the modal response quantities |q| <|q|>, lal...» and 


max’ 
<|q|> are obtained and multiplied times the appropriate modal hydrodynamic 
</Q|>, [Q| 


the problem arises of how to combine the modal responses. 


parameter participation factor to obtain [Q| , and </Q|>, 


max? max 


To obtain an estimate of {|x } and {|x }, the modal responses 


aay 5 hee 


can be combined in any one of several methods as is the case for obtaining 


{|Z| 


 - An upper bound is obtained by summing the absolute values of the 
modal responses (ZABS). 

An upper bound on the quantities {<|x.|>} and {<|x_|>} can also be 
obtained by the ZXABS method of modal combination, but this is not readily 
apparent. That the XABS method of modal combination yields an upper bound 


on the quantity {<|x.|>} (or {<|x,|>}) is shown as follows: 
{x} = [o]{Q} (3.42) 


Assume that the vector {Q(t)} has been obtained, then, at the jth node of 


the system, 


x (8) = 651 Q Ct) # 652%(E) + «05 amt (3.53) 
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Taking the absolute value of both sides of Eq. 3.53 


Q(t) | (3.54) 


jJsnm‘nm 


(t)] = [45,9)(4) + o,,Q(t) +... @ 


An upper bound on IX. (t)| is seen from Eq. 3.54 to be 
J 


aad Sos LIQ Q(t] + Loso]1Qo(t)] + ~ © © 1% ml 1G pm t | 
(3.55) 
Integrating both sides of Eq. 3.55 over the time of response, 
ae Tea 
t | Ix, (lat s £ ] Chop llQy(edl + LojghiQct)| 
0 0) 
Pe a Oral hood dt (3.56) 
Or, 
six, I < los 1<1Q) 1 ni Ido l<1Qo1> 
Da MES aed) Oa (3.57) 


Which is seen to be the XABS method of modal combination. 
Rules for constructing idealized hydrodynamic parameter response 


spectra from which |q| , and <|q|> may be estimated for each 


max’ <|q|>, [4 


<]Q|>, |Q| 


the hydrodynamic parameter participation factors) of the multi-degree of 


max 


mode (and then |Q| , and <|Q|> are computed from the q's and 


max? max 


freedom system are discussed in the next section. 
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3.5 Idealized Elastic Hydrodynamic Parameter Response Spectra 


Idealized elastic hydrodynamic parameter response spectra can be 


constructed for each mode of the multi-degree of freedom system in the fol- 


lowing manner: 


1. 


Compute the natural frequencies, ['w-], and mode shapes, [4], 
for the system. 
Establish the structural damping ratio, B. > for each mode. 


1 
Compute By and HDPF . from 
| 


nn 
o2\° 
_ | 97] Dj aia 22 
_ nn ' Bee. (oe ) 
ahead 
[28,0] = acy pele] (3.24) 
{HDPF} = peace (3.49) 


With the quantities Waa Bs 5° Baa? Mm (or <[XQ|>) 


X 
g max > “g max 


(or <|x0|>), and x (or SNkele) construct the idealized 


g max 
elastic hydrodynamic parameter response spectrum for the ith 
mode exactly as described in Section 2.7 for the single degree 


of freedom system. 


When an idealized spectrum has been drawn for each mode of the sys- 


tem, maximum modal hydrodynamic parameters are obtained from the spectra, 


multiplied times the appropriate modal hydrodynamic parameter participation 


factor, and combined in some manner to yield an approximate solution of the 


hydrodynamic parameter of interest. 
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CHAPTER 4 
APPLICATIONS 
4.1 General 


In Chapters 2 and 3 a method was developed for linearizing the non- 
linear hydrodynamic drag force on a single degree of freedom oscillator and 
a multi-degree of freedom lumped parameter model. The method has been shown 
to yield satisfactory results for systems subjected to earthquake excitation 
of the type used in this thesis. Although the method of linearizing employed 
in the multi-degree of freedom system places a linear damper at the level of 
the top mass where a nonlinear damper does not exist, it has been shown to 
yield maximum response quantities which are within a very few percent of the 
response quantities computed for the nonlinear system. Further, linearizing 
the drag force by using velocity characteristics of the ground motion has been 
Shown to yield generally conservative results in all but the low frequency 
range. The linearization technique has been developed in such a manner that 
the resulting linear simultaneous differential equations of motion for the 
multi-degree of freedom system possess classical normal modes of response. 
The existence of classical normal modes permits the use of normal mode 
Superposition, with the modal response quantities determined by response 
spectrum techniques, aS an approximate method of computing the response of 
multi-degree of freedom systems. The approximate method can be used for dis- 
placement response quantities and hydrodynamic parameter response quantities. 
The research reported in this thesis is not intended to reinvent the method of 


solution by the superposition of normal mode solutions as determined by 
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response spectra, but is devoted primarily to developing the linearizing 
technique which insures the existence of classical normal modes and yields 
generally conservative results. 

This chapter demonstrates the use of the principles and methods 
developed and discussed in Chapters 2 and 3 in computing the response of 
a series of actual offshore tower designs and one of the test systems des- 


cribed in Chapter 3. 
4.2 Displacement Response 


Table 4.1 lists the properties of a series of offshore towers de- 
signed for water depths of 400, 600, 800, and 1000 feet.* A structural damp- 
ing matrix [C.] which yields modal damping ratios of .01 in each mode is com- 
puted for each tower as described in Chapter 3. The structures are subjected 
to the El Centro earthquake depicted in Fig. 1.2, and various response quan- 
tities are computed. In the first set of computations, the equations of mo- 
tion are numerically integrated and response quantities are found for the 
case in which there are no hydrodynamic drag forces acting on the structure. 
The same response quantities are next computed for the case in which the actual 
nonlinear drag forces are acting on the structure. Finally, the nonlinear 
drag forces are linearized by using the velocity characteristics of the 
ground motion, and the response quantities are computed by the combination 
(ZABS and SRSS) of maximum modal responses found from the true response spec- 
trum for each of the first five modes of the structures. The amount of linear 
hydraulic damping in each mode found from the linearization procedure is shown 


in Table 4.2. The response quantities of interest in these calculations are 
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the maximum displacement between the top deck and the ground, the maximum over- 
turning moment, and the maximum base shear expressed as a percent of the total 
tower weight, including the added mass of the structure. The results of the 
calculations are shown in Fig. 4.1. 

A study of Fig. 4.1 reveals that, for this series of towers sub- 
jected to the El Centro earthquake, the nonlinear hydrodynamic drag forces 
reduce the response of the tallest tower in the top deck displacement and 
overturning moment quantities by approximately 30 percent from the "no drag" 
case. The maximum base shear for the tallest tower is reduced only about 8 
percent by the nonlinear hydrodynamic drag forces. Exact response calculations 
were made for the other three earthquake records used in this thesis and the 
reduction in response due to drag was found to be somewhat lower in magnitude 
than for the EI Centro earthquake. For the tower in 1000 feet of water, the 
reduction (due to hydrodynamic drag) in the maximum top deck displacement 
relative to the ground was found to be about 8.5 percent for the Taft earth- 
quake, 11.5 percent for the Ferndale earthquake, and 17 percent for the San 
Fernando earthquake. In general, the taller, more flexible towers showed 
more reduction in response due to the hydraulic drag than did the shorter, 
stiffer towers. No generalizations can be made, however, regarding the mag- 
nitude of the reduction in response relative to tower height or stiffness 
Since the actual response is such a complex combination of variables in both 
the structure and the earthquake. 

In the superposition of maximum modal responses computed from the 
true response spectrum for each of the first five modes of vibration, it was 


found that the higher modes contributed rather strongly to the maximum base 
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shear. The contribution of each mode to the maximum base shear computed by 
the ZABS combination of modal responses is shown in Table 4.3. Although the 
SRSS combination of modal responses appears from Fig. 4.1 to provide a very 
good estimate of the response of this series of towers to the El Centro earth- 
quake, generalities concerning the accuracy of one method of modal combination 
relative to another method of modal combination cannot be made.’? 

The response of the towers is computed by the combination of maxi- 
mum modal responses (ZABS and SRSS) found from the idealized deformation re- 
sponse spectrum for each of the first five modes of vibration. The idealized 
spectra were constructed in accordance with the procedure described in Chapter 
3 for ground motion quantities equal to those of the El Centro earthquake. The 
idealized spectra for the first five modes of tower C are shown in Fig. 4.2. 
The estimated hydraulic damping ratios in the 3rd, 4th, and 5th modes of the 
structure are such that the resulting idealized spectra for those modes cannot 
be shown separately and distinctly on the scale of Fig. 4.2. It is noted that, 
for the idealized spectra shown in Fig. 4.2, the frequencies of the system are 
well removed from the high frequency tail region of their respective response 
spectrum. This was found to be the case for each of the towers A, B, C, and 
D. The results of the computations using the idealized spectra are shown in 


mg. 4.3. 
4.3 Hydrodynamic Parameters 


To demonstrate the computation of the hydrodynamic parameters 


ae <[x.|>; Xela, and <|x.|> at each node of a multi-degree of freedom 


System, the five degree of freedom system 4M is used. The properties of this 
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system, with a structural damping matrix that yields a damping ratio of .02 
in each mode, are shown in Table 4.4. The magnitude of [K] was chosen to 
yield frequencies of the system that fall into the velocity amplification 
region of the spectrum. The hydrodynamic parameter response quantities are 
computed by the numerical integration of the equations of motion and by the 
combination (ZABS and SRSS) of maximum modal responses found from the true 
spectral response in each mode. The results of the computations are shown 
in Table 4.5. It is seen from Table 4.5 that the combination of modal re- 
sponses provides results which are, in general, sufficiently close to the 
exact response at each node to enable the selection of inertia and drag 
coefficients C and Cy with reasonable confidence. In lieu of using the 
true spectral response quantities, which require the use of a digital com- 


puter to obtain for each mode, an idealized spectrum can be constructed for 


each mode in the manner set forth in Chapter 3. 
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CHAPTER 5 
CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 
5.1 Conclusions 


In this thesis, an approximate method has been developed for com- 
puting the response of a multi-degree of freedom system, partially submerged 
in water, to earthquake excitation. In Chapter 2, the effects of the non- 
linear hydrodynamic drag force on a single degree of freedom oscillator are 
shown, and a generally conservative method of linearizing the drag force is 
developed. Amplification factors are developed which permit the construction 
of idealized response spectra for deformation and hydrodynamic parameters. 
Chapter 3 extends the principles from the single degree of freedom oscillator 
to multi-degree of freedom systems and develops a generally conservative 
method of linearizing the drag forces. The method of linearizing insures the 
existence of classical normal modes in the equations of motion, which in turn 
permits the response quantities of interest to be obtained by the superposition 
of normal mode responses. In Chapter 4 the response to the El Centro earthquake 
of a series of offshore towers is computed using the methods developed in 
Chapter 3. 
Several conclusions may be drawn from the development of principles 
and methods and the results of sample problems in the preceding three chapters. 
1. The maximum relative displacement response of a single degree 
of freedom oscillator, with drag characteristics described by 
KD ‘ max 
K 


the parameter » approaches a constant value at high 


frequencies as a result of the nonlinear hydrodynamic drag 
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force growing very large relative to the inertia force in the 
equation of motion. The magnitude of this relatively constant 
high frequency tail can be estimated for earthquake motions by 


the expression 


- °2 
Z| 7 KD Xg max xg max 
max — K Xq ae 











Linearizing the drag force on a single degree of freedom oscil- 
lator by using the absolute velocity response of the oscillator 
yields slightly conservative results for the maximum response 

to earthquake motions until the high frequency, high drag tail 
begins. 

Linearizing the drag force on a single degree of freedom oscil- 
lator by using the velocity characteristics of the ground motion 
yields conservative results for the maximum response to earth- 
quake motions except for low frequency systems and very high 
frequency systems. 

The conservatism introduced by using the ground motion velocity 
characteristics in the linearization of the drag force is re- 
duced as the amount of structural damping in the system is in- 
creased. 

The maximum response of the single degree of freedom oscillator 
to earthquake motions is sensitive to the total amount of linear, 
velocity dependent damping in the system, and not to the distri- 
bution of this damping between structural and hydraulic, up to 


linear hydraulic damping ratios of 10 percent of critical. 





ae 
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The nonlinear hydrodynamic drag has very little effect on the 
elasto-plastic deformation response spectrum in the low fre- 
quency region. 

A high frequency tail exists in the elasto-plastic deformation 
response spectrum for a system with nonlinear hydrodynamic drag 
forces. Large amounts of drag cause greater reductions in the 
yield deformation, zy from the maximum elastic deformation for 
a given ductility factor than do small amounts of drag. 

The nonlinear hydrodynamic drag force has a damping effect on 
the elasto-plastic deformation response spectrum in the mid- 
frequency range. The damping effect of the drag becomes less 
pronounced as the ductility factor increases. 

Idealized response spectra which consider the effects of fluid- 
structure interaction can be constructed as described in this 
thesis. 

The linearization of the drag forces on a multi-degree of freedom 
system by using the absolute velocity response at each node and 
forcing the linear drag factor matrix [Ky ~J to be proportional 
to the mass matrix [-M.] as described herein is satisfactory. 
Linearizing the drag forces on a multi-degree of freedom system 
by using the velocity characteristics of the ground motion and 
forcing the drag factor matrix to be proportional to the mass 
matrix as described herein yields generally conservative results 
except for low frequency systems. 


For the multi-degree of freedom systems studied in this thesis, 
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14. 


Foe 
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the response to earthquake excitation is relatively insensitive 
to the spatial distribution of the linear drag forces. 

The hydrodynamic drag forces can cause a significant reduction 
in the response of lightly damped systems to earthquake motions. 
The amount of reduction in response due to hydrodynamic drag 
varies considerably from one response quantity to another and 
from one earthquake to another. 

The reduction in response due to hydrodynamic drag is decreased 
as the amount of structural damping in the system is increased. 
Higher modes of a multi-degree of freedom system can contribute 


appreciably to the total response of the system. 


5.2 Recommendations for Further Study 


The approximate methods of analysis developed in this thesis are 


limited by the types of multi-degree of freedom systems used in the investi- 


gation and by the inherent limitations in the hydrodynamic force equation 


used. Further, the study in this thesis is devoted to excitation of the 


structure in still water, and the ocean is never completely at rest. There 


are several areas in which further studies should be made of the response of 


offshore structures to earthquake excitation. 


1. 


The validity of the lumped parameter model should be tested 
through the use of field measurements of the response of actual 
offshore structures to strong motion earthquakes. 

Model tests of single degree of freedom oscillators and multi- 


degree of freedom structures should be conducted for the types 
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of base motion considered in this thesis. Such model tests will 
be useful in verifying the principles developed herein. 

Since the ocean is never actually at rest, analytical and ex- 
perimental studies should be conducted on the response of sub- 
merged and partially submerged structures subjected to earthquake 
motions combined with forces resulting from waves and currents. 
The studies in this thesis assume ground motion in one direction 
only. It is recommended that studies be made of the response 
when the two orthogonal components of the horizontal ground motion 
are acting on the structure at the same time. 

The earthquake response of the structure in torsion and in the 
vertical direction should be investigated. 

It is recommended that the effects of foundation and soil con- 
ditions be included in the model. There are several reasonable 
approaches to modeling the foundation which should be investi- 
@ated,-° °° 

The hydrodynamic force equation used in this thesis (Morison's 
equation) is based on studies with relatively rigid piles. There 
is some evidence’?>*° that the forces on flexible, oscillating 
Structures may be better expressed through the use of different 
drag and inertia coefficients, or even through the use of a 
different force expression. Much more information is needed 
concerning the hydrodynamic forces acting on flexible, oscillating 
structures. 


The linearizing technique should be tested for suitability of 
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use when the structure has radically different distribution of 
mass, stiffness, and drag properties from those test structures 
used in this thesis. 

It is recommended that the effects of the hydrodynamic forces 

on the inelastic response of the structure be studied in greater 
depth. While design for elastic response is satisfactory for 
earthquakes which may be expected to occur frequently, elastic 
design for infrequent, severe earthquakes may result in ex- 
cessively conservative designs. Some response into the in- 
elastic range should be permitted for the occurrence of an un- 
usually strong earthquake, and approximate methods are desirable 
for the estimation of this inelastic response. 

It is recommended that approximate methods be developed for 

the estimation of the response to stochastic representations 


of the forces resulting from earthquakes, waves, and currents. 
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Table 2.1] 


Characteristics of Series Representation 
of Base Acceleration 


Quantity Maximum Amplitude Time of Occurrence 
Xa Z56c6 0235 ty 
Xa 3.65 Tt€/ Sec 0.42 ty 
Xq OeZ071 0.50 ty 


xX. = -0.20401 cos 2mt - 6.87108 cos 4mt + 38.35016 cos 6nt 
= dro 5225 GOS ont = croucoec cos 10nt 
+ 12.91937 cos l2nt 


(X in ft/sec?) 
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Table 2.3 


Distribution of Total Damping Between Structural 
and Hydraulic for Comparison Purposes 





Bh For 
total Be = 0 B. = .02 B. = .05 B. = .10 
.02 .02 0 
.05 “05 0 
.07 705 .02 
a0 mae) .05 0 
BZ a0 .02 
aS oO AOS 
Table 2.4 


Elastic Deformation Spectrum Amplification Factors 


Amplification Factor 


Be otal \*) Displacement Velocity Acceleration 
0 eS 4.0 6.4 
OFS ae 520 556 
| 2.0 Sys ssa 
2 1.8 2.8 4.3 
5 1.4 i238 2.6 
/ lec Vets 1.9 
10 1.1 tas L25 
20 1.0 ia eZ 
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Table 2.5 


Elastic Hydrodynamic Parameter Spectrum 
Amplification Factors 


paneer Spectrum 
Amplification Factors 
Ptotal Displacement Velocity 
(%) 
0 2.4 52 
] Cae 4.4 
2 Zeal Sia 18) 
5 1.8 Sell 
10 le Za 
<|x.|> Spectrum 
Amplification Factors 
PLotal Displacement Velocity 
(%) 
0 Sine: 10.9 
] 3.4 Sa 
a oO 6.8 
5 ) 33 4.2 


10 ere: Sad 





Btotal 


(%) 


Brotal 
(%) 


Table 2 


.5 (Continued) 


al ie Spectrum 


Amplification Factors 


Displacement 


Lae 


1.4 
WZ 


Displacement 


3.4 
on) 
2./ 
Zaedl 
ete 


Velocity Acceleration 
4.1 fea! 
gol 4.9 
2.6 oo 
e720 cee 
1.8 2.6 


<|x.|> Spectrum 
Amplification Factors 


Velocity Acceleration 
10-7) 13.5 
6.3 6.2 
5] a8 
3.0 2.9 


cay eae 
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Table 2.6 


Actual and Approximate Reduction Factors for the 
High Frequency Region of Inelastic Response 
Spectra. E] Centro Earthquake 


KX 
B. = 0.02, a 0.50 





Actual 

Ductility Factor Z, felts oe mF = 
u (inches) (inches ) y 

I 0.63676 0.63676 1.00 1.00 

168 0.49526 0.74289 On77 Dev 

2 0.42835 0.85670 0.67 0.67 

3 0.35975 1.07925 0.56 0.60 

4 Des 1523 1.26093 0.49 257 


Note: All calculations are at a frequency of 10 cps. 
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Estimated Linear Hydraulic Damping Ratios 
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Fig. 2.24 Elastic Deformation Response Spectra for San Fernando Earthquake 
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